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I.  Introduction 

The  objective  of  the  work  described  in  this  report  is  to  use  the  NSWC  Mathematical 
Library  (MLIB)[1],  which  is  written  in  FORTRAN  77,  as  an  aid  in  testing  a  new  Fortran 
90  compiler  used  with  the  CRAY  EL98  mainframe  computer.  Indirectly,  one  can  view  the 
testing  as  an  aid  in  establishing  whether  MLIB  can  operate  in  a  Fortran  90  environment. 
The  final  result  is  a  transportable  driver  program,  TEST90,  which  is  not  restricted  to  testing 
only  the  new  90-compiler.  Any  post  Fortran  77  compiler  can  be  tested  provided  it  is  limited 
to  a  single  precision  word  length  of  64  bits  or  less.  This  limitation  is  imposed  because  the 
routines  in  MLIB  are  based  on  a  64  bit  word  length  with  a  48  bit  floating  point  mantissa. 

TEST90  has  been  used  to  test  the  latest  Fortran  90  compilers  of  the  CRAY  EL98  and 
IBM  PC  machines;  no  errors  were  found,  however  TEST90  did  reveal  a  complex  arthmetic 
error  in  an  earlier  Cray  EL98  compiler.  MLIB  routines  ran  under  TEST90  without  any 
problems  on  both  machines. 

MLIB  contains  1099  routines  of  which  586  are  callable  directly  by  the  user.  The  rest 
are  called  supporting  routines.  We  have  chosen  57  of  the  callable  routines,  Group  A,  as  a 
representative  sample  of  MLIB.  Eleven  other  callable  routines  are  included  in  Group  A  that 
act  as  supporting  routines  making  a  total  of  68.  In  addition,  17  double  precision  versions 
of  a  subset  of  the  57,  which  are  not  part  of  MLIB,  were  constructed  using  a  routine  called 
FETCH  (FCH),  (designed  by  Russ  Gnoffo.  See  CD  disk  of  [1]).  These  routines  are  said  to 
belong  to  Group  B.  Thus  85  routines,  making  up  Groups  A  and  B,  were  used  for  testing. 
Supporting  these  routines  were  approximately  282  additional  routines. 

Thirteen  hundred  and  twenty-five  output  quantities  were  generated  and  examined.  The 
fixed  input  data  specified  for  these  test  cases  were  chosen  to  bring  into  play  as  many  of  the 
supporting  routines  as  possible.  Of  the  282  supporting  routines  mentioned  in  the  previous 
paragraph,  238  (or  319  of  the  total  367)  were  actually  executed  in  running  the  1325  test 
cases. 

TEST90  also  calls  subprograms  which  check  that  all  the  bit  and  string  commands  of  the 
Fortran  90  compiler  operate  correctly.1 

For  the  numerical  testing  of  groups  A  and  B,  a  true  result  for  each  output  was  obtained 
independently  to  35  significant  digits  by  using  Maple V  software  [2].  TEST90  checks  that 
each  Fortran  90  test  case  calculated  result  agrees  with  the  corresponding  true  Maple  value 
to  within  a  certain  accuracy.  In  order  to  do  this,  for  each  single  precision  (double  precision) 
output,  a  prespecified  tolerance  is  given  in  terms  of  a  positive  parameter  EPS  (EPD)  that 

1MLIB  does  not  use  any  bit,  string,  or  I/O  statements. 
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depends  on  the  machine  dependent  parameter  SPMPAR(l)  (DPMPAR(1))[1,  p.3-4]2.  Then 
if  the  relative  error,  which  is  defined  in  (1),  is  within  the  tolerance,  that  case  is  said  to 
have  passed.  If  the  relative  error  exceeds  the  tolerance,  a  message  is  printed  by  TEST90 
identifying  the  routine,  the  case  number,  and  the  output;  that  particular  test  case  has  failed. 
If  no  disagreement  occurs  for  all  1325  cases,  a  message  is  printed  at  the  end  of  the  run  that  all 
has  gone  well.  It  should  be  noted  that  in  case  of  a  failure  the  test  case  and  the  subprogram 
of  Group  A  or  B  is  identified,  but  the  specific  problem  with  the  compiler  is  not  isolated.  An 
independent  study  is  required  to  determine  what  caused  the  failure. 

Let  T  denote  the  true  result  (obtained  from  Maple),  C  the  calculated  result  (one  of  the 
1325  numerical  outputs),  and  let,  for  asingle  precision  result,  X  =  SPMPAR(2)/SPMPAR(1) 
and  for  double  precision,  X  =  DPMPAR(2)/DPMPAR(1).  Then  the  relative  error  is  given 
by 

r  max(|T|,  |C|)  for  |T|  <  X. 

RE  =  l  |T|  for  |C|  <  X.  (1) 

[  Q  for  Q  >  5  SPMPAR(1)(DPMPAR(1))  else  0, 

where  Q=  |(T  —  C)/T|. 

A  test  result  is  said  to  pass  if 

RE  <  EPS  (or  EPD  for  a  double  precision  result),  (2) 

where  the  machine  dependent  parameters  EPS  and  EPD  are  defined  below. 

It  is  convenient  here  to  introduce  the  following  briefer  notation: 


CRAY: 

SPMPAR(l) 

— ♦  Sc, 

DPMPAR(l) 

Dc, 

RE 

REc 

IBM  PC: 

SPMPAR(l) 

-+  Spc, 

DPMPAR(l)  -> 

Dpc, 

RE  -> 

REpc 

Machine: 

SPMPAR(l) 

— >  Sm, 

DPMPAR(l)  -*■ 

Dm, 

RE  -4 

REm 

where  Machine 

refers  to  the  computer  which  is  used  with  the 

compiler  being  tested. 

The  quantities  EPS  and  EPD  are  defined  in  terms  of  pre-computed  elements  SE  and  SF, 
where  SE  is  given  by 

{REc/Sc,  if  REc  >  5  Sc;  else  1.25  (single  precision). 

REc/Dc,  if  REc  >  5  Dc;  else  1.25  (double  precision).  (3) 

1,  if  T  =  C  =  0. 

2SPMPAR(1)  (DPMPAR(l))  denotes  the  smallest  positive  value  such  that  1  +  SPMPAR(DPMPAR(1))  >  1. 

CRAY  :  SPMPAR(l)  (DPMPAR(l))  =  2~*7  (2~95).  PC  :  SPMPAR(l)  (DPMPAR(l))  =  2"23  (2"62). 

3SPMPAR(2)  (DPMPAR(2)),  is  the  smallest  positive  single  precision  (double  precision)  number  in  the  machine. 

CRAY  :  SPMPAR(2)  =  DPMPAR(2)  =  2~8188;  PC  :  SPMPAR(2)  (DPMPAR(2))  =  2~126  (2~1021). 
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Note  that  the  1325  SE  values  depend  only  on  the  64  bit  CRAY. 

In  order  to  define  SF  (safety  factor)  for  a  particular  output,  it  is  necessary  to  determine 
whether  rounding  error  TZ  or  truncation  error  T  is  the  major  contributor  to  the  relative  error. 
For  example,  in  the  inversion  of  the  Hilbert  matrix  TZ  plays  the  dominate  role,  whereas  in  an 
application  of  Runge  Kutta,  T  is  the  main  contributor.  With  this  in  mind,  each  of  the  1325 
cases  is  classified  as  either  an  TZ  or  a  T  type.  Note  that  this  classification  is  independent  of 
the  machine.  Thus  let 


SE  *  Spc(Dpc), 
SE  *  Sc(Dc), 


if  REpc  is  of  type  TZ  (single  (double)  precision) 
if  REpc  is  of  type  T  (single  (double)  precision). 


Using  the  above  parameters,  we  have 


{5,  if  REpc  <  5  Spc(Dpc)  (single (double)  precision) 

Max[3, 5  REpc/Z  ],  otherwise. 


(4) 

(5) 


Hence  SF  is  computed  using  the  32-bit  PC,  as  well  as  SE  which  was  computed  using  the 
CRAY.  This  approach  allows  TEST90  to  be  used  on  computers  with  word  lengths  from  32 
to  64  bits. 

Using  SE  and  SF,  from  (3)  and  (5),  the  tolerance  parameters  EPS  and  EPD  are  given  by 


EPS  = 


Max[ 5  *  Sm,  SF  *  SE  *  Sm  ],  if  REpc  is  of  type  TZ 
Max[5  *  Sm,  SF  *  SE  *  Sc],  if  REpc  is  of  type  T. 


(6) 


EPD  = 


Max[  5  *  Dm,  SF  *  SE  *  Dm  ],  if  REpc  is  of  type  TZ 
Max[5  *  Dm,  SF  *  SE  *  Dc],  if  REpc  is  of  type  T. 


(7) 


If  a  single  precision  output  result  is  tested  using  (2),  where  EPS  is  given  by  (6),  it  passes  if 
REm<EPS  or  if  EPS  >  1/10  for  TZ  type  cases.  The  latter  inequality  holds  when  no  digits  in 
C  are  correct  because  of  using  too  short  a  word  length  (An  example  is  attempting  to  find  the 
inverse  elements  of  the  tenth  order  Hilbert  matrix  in  single  precision  on  a  32  bit  machine). 
If  a  double  precision  output  is  tested,  using  (7),  it  passes  only  if  REm<EPD.  All  1325  cases 
will  pass  if  no  compiler  errors  are  found. 

Section  II  lists  the  call  lines  of  the  Group  A  and  the  Group  B  routines,  the  associated  test 
cases,  and  the  corresponding  35  digit  true  results  obtained  from  Maple  software.  Subroutines 
testing  bit  manipulation  functions  are  described  in  Section  III  and  subroutines  testing  string 
operations  are  explained  in  Section  IV.  Section  V  contains  directions  for  using  TEST90. 
Section  VI  contains  an  index  of  the  routines  discussed  in  Section  II. 
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II.  Callable  Routines  and  True  Value  Maple  Data 

In  this  section  we  list  the  directly  callable  routines  that  were  used,  with  their  call  line 
from  MLIB  and  their  page  number  in  [1],  and  the  routines  in  Group  B.  Also  given  are  the 
test  cases,  and  the  35  digit  Maple  [2]  results  which  are  considered  to  be  the  true  values. 
Each  of  the  subprograms  is  numbered  in  numerical  order,  1  <  N  <  74,  and  on  the  same  line 
in  parentheses  the  range  of  each  set  of  outputs  is  itemized  in  numerical  order,  so  that  on 
the  first  line  with  each  subprogram  we  have  N.  (ajj  —  b^)  —  PageM,  where  1  <  a^  <  bu  < 
1325,  aN+i  =  bN  +  1,  and  M  refers  to  the  page  number  in  [1].  A  complex  number  is  listed 
as  an  ordered  pair  of  real  numbers,  i.e.,  z  =  a  +  lb  =  (a,  b);  I  =  yf—  1.  If  the  imaginary  part 
of  a  complex  number  is  not  given  it  is  assumed  to  take  the  value  zero. 
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1.  (1-7) -Page  11 

y  =  CBRT(x)  -  CUBE  ROOT. 

xi,yi  =  1E(— 200), 

.21544346900318837217592935665193505E(— 66) 

x2,  y2  =  5E(— 29), 

.36840314986403866057798228335798072E(— 9) 

x3,y3  =  1E(— 3), 

.1 

*4,  y4  =  0, 

0 

x5,  y5  =  1E3, 

10 

xg,  y6  =  5E29, 

.7937005259840997373758528196361 541 3E1 0 

x7,  y7  =  1E50, 

.4641 58883361277889241 007635091 94466E1 7 

2.  (8— 14)— Page  11 

y  =  DCBRT(x)  -  DOUBLE  PRECISION  CUBE  ROOT. 

Double  precision  version  of  CBRT.  Uses  same  input. 

3.  (15— 18)— Page  11 

w  =  ARTNQ(y,  x)  -  FOUR  QUADRANT  ARCTANGENT,  [0,  2tt). 
yi/xi  =  3.9/1. 7,  Wl  =  1.1597317794050322209919342332641325 

y2/x2  =  3.9/ -1.7,  w2  =  1.9818608741847610174707091500153704 

y3/x3,  =  -3.9/  -1.7,  w3  =  4.3013244329948254594545776165436354 

y4/x4,  =  — 3.9/1.7,  w4  =  5.1234535277745542559333525332948733 

4.  (19— 22)— Page  11 

w  =  DARTNQ(y/x)  -  FOUR  QUADRANT  ARCTANGENT. 

Double  precision  version  of  ARTNQ.  Uses  same  inputs  as  ARTNQ. 

5.  (23— 30)— Page  13 

Call  DCSQRT(Z,  W)  -  n/Z,  DOUBLE  PRECISION  COMPLEX. 

Z  =  (real,  imag),  Za  =  (0, 1),  Z2=(l,l),  Z3  =  (50, 25),  Z4  =  (1E30, 1E40). 

W1=(.70710678118654752440084436210484905,  .70710678118654752440084436210484905) 
W2=  (1 .0986841 1 34678099660398011952406784,  .45508986056222734130435775782246858) 
W3=  (7.27673345 1 1 267740406 1 3309 1 985484853 ,  1.7178037486125623206928287195727927) 
W4=(70710678122190286346.105562180133025,  70710678115119218534.240086936133420) 
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6.  (31 -35) -Page  29 

y  =  REXP(x)  -  [EXP(x)  -  1]. 

Xl  =  -50,  x2  =  -5,  x3  =  1E(— 20),  X4  =  50,  x5  =  85. 
yi  =  -.99999999999999999999980712501520361 

y2  =  -.99326205300091453290336395157685158 

y3  =  .10000000000000000000050000000000000E(— 19) 

y4  =  .51847055285870724640864533229334854E22 

y5  =  .822301 27146229135 1 030432801 64077747E37 

7.  (36— 40)— Page  29 

y  =  DREXP(x)  -  [EXP(x)  -  lj. 

Double  precision  version  of  REXP.  Uses  same  input. 

8.  (41-46)— Page  31 

y  =  RLOG(x)  -  [x  -  1  -  LOG(x)]. 

Xi  =  [lE(-lOO),  1E(— 5),  1,  1.005,  50,  100];  i=l,...  ,6. 

yi  =  229.25850929940456840179914546843642 

y2  =  10.512935464970228420089957273421821,  y3= 

y4  =  .12458488960926387897797540656528063E(— 4) 

y5  =  45.087976994571853941381249212089448 

y6  =  94.394829814011908631964017090631272 

9.  (47-52) -Page  32 

yi  =  DRLOG(xj)  -  [x  -  1  -  LOG(x)]. 

Double  precision  version  of  RLOG.  Uses  same  input. 

10.  (53— 59)— Page  46 

y  =  ERF(x)  -  ERROR  FUNCTION. 

x3=  [-5,  -1,  .4,  1.3,  3.5,  6.7,  20];  i  =  l,...  ,7. 
yi  _  -.99999999999846254020557196514981166 

y2  =  -.84270079294971486934122063508260926 

y3  =  .42839235504666845510360384532017244 

y4  =  .93400794494065243660389332750371102 

y5  =  .99999925690162765858725447631624390 

y6  =  .99999999999999999999733828542362341 

yr  =  1-0 
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11.  (60-66) -Page  48 

y  =  DERF(x)  -  DOUBLE  PRECISION  ERROR  FUNCTION. 

Uses  same  arguments  as  ERF. 

12.  (67-73) -Page  46 

y  =  ERFCl(IND,x)  -  ex2[l  -  ERF(x)]  =  ex2ERFC(x). 

Uses  same  arguments  as  ERF.  IND  =  1. 
yx  =  .14400979867466104041058963430588211E12 

y2  =  .50089800807622834663098245982148099E1 
y3  =  .67078778529476152332941247446792512 

y4  =  .35764266908609031764697338728600749 

y5  =  .15529365560889429740272649758187893 
y6  =  .832991 894668104129940674394787281 64E(— 1 ) 

y7  =  .281 74348741 051 31 931 86491 54534470758E(-1) 

13.  (74-80) -Page  49 

y  =  DERFC1(IND,  x)  -  DOUBLE  PRECISION  ERFC1. 

Uses  same  arguments  as  ERF.  IND  =  1. 

14.  (81— 88)— Page  45 

Call  CERF(MO,z,w) -ERROR  FUNCTION,  COMPLEX  ARGUMENT. 

Zi  =  [(.4,  .3),  (3.4,  2.3),  (5.3,  3.7),  (6.2,  .009)];  i  =  1, . . .  ,  4;  MO  =  0. 

wi  =  (.46443721297956140209712031572604685,  .29443980776931707303974297404552757) 

w2  =  (1.0002253313625125573357082824478311,  -.12334728086728759257490539663821289E(-3)) 

w3  =  (1.0000000252590040600802224870937450,  .41332248953387653013128019073968635E(-7)) 

w4  =  (.99999999999999999819476945857815225,  ,20489319058601197920097004639286141E(-18)) 

15.  (89— 96)— Page  47 

Call  DCERF(MO,  Z,  W)— CERF,  DOUBLE  PRECISION  COMPLEX. 

Uses  same  arguments  as  CERF.  MO=0. 

16.  (97-1 04) -Page  46 

Call  CERFC(MO,z,  w)  -  ez2ERFC(z)  WITH  COMPLEX  ARGUMENT. 

Uses  same  arguments  as  CERF.  MO=l,  Re(z)>  0. 

wi  =  (.63299603234343978026485063394958694,  -.17020263553343032296655297457404080) 

w2  =  (.11438075139307192540234735774672486,  -,73122766392201685536727808849318310E(-1)) 

w3  =  (.71 8000861 51 584466507931 675846826108E(— 1 ) ,  -,48959427476482652124106743611245310E(-1)) 
w4  =  (.89857930230247263981968919036797811E(— 1),  -.1272473689470905321 9621 092389506531E( -3)) 
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17.  (105-1 12)-Page  48 

Call  DCERFC(MO,  Z,  W)— CERFC(z),  DOUBLE  PRECISION. 
Uses  same  arguments  as  CERF.  MO=l,  Re(Z)>  0. 


18.  (113-184)— Page  88 

Call  BRATIO(a,  b,x,y,w,wl,  IERR) -INCOMPLETE  BETA  FUNCTION  . 
x  +  y  =  1,  Output :  w,  wl  (w  +  wl  =  1),  IERR(=  0), 

A=ai=[.25, 1.25,  7.25,  20.25],  B=bj=[.25,  2.25,  8.25,  25.25],  X=xk=[.25,  .5,  .75,  .99] 

All  combinations  of  the  elements  of  A,B,X  are  used.  Data  below  shown  as  BRAT(A[i],  B[j], 
X[k])=  w,  or  BRAT(B[j],  A[i], y[k])=wl, where  the  first  is  listed  if  w<  wl  otherwise  the 
second  is  listed.  For  i=l,  j=2,  k=3,  Call  BRATIO(.25,  2.25,  .75,  w,  wl,  IERR)  so 
that  BRAT(2.25,  .25,  .25)  =  wl  =  .7351410084060216157635651 1063802129E(-2).  In 
addition,  another  8  values  of  w  =  BRAT(a,  b,  x)  are  included  at  the  end  of  the  listing 
that  follows,  making  a  total  of  72  BRAT  cases. 


BRAT(A[1],  B[l],  X[l]) 
BRAT(A[1],  B[l],  X[2]) 
BRAT(B[1],A[1],1-X[3]) 
BRAT(B[1],  A[l],  1-X[4]) 
BRAT(B[2],  A[l],  1-X[1]) 
BRAT(B[2],A[1],1-X[2]) 
BRAT(B[2],  A[l],  1-X[3]) 
BRAT(B[2],  A[l],  1-X[4]) 
BRAT(B[3],  A[l],  1-X[1]) 
BRAT(B[3],A[1],1-X[2]) 
BRAT(B[3],  A[l],  1-X[3]) 
BRAT(B[3],  A[l],  1-X[4]) 
BRAT(B[4],A[1],1-X[1]) 
BRAT(B[4],A[1],1-X[2]) 
BRAT(B[4],A[1],1-X[3]) 
BRAT(B[4],A[1],1-X[4]) 
BRAT(A[2],  B[1],X[1]) 
BRAT(A[2],B[1],X[2]) 
BRAT(A[2],B[1],X[3]) 
BRAT(B[1],  A  [2],  1-X[4]) 


=  .39775677831735583689234908829962206 

=  .50000000000000000000000000000000000 

=  .39775677831735583689234908829962206 

=  .17081535717490444338991136891551120 

=  .14085505922511297496662594047200120 

=  .4234414189891208716260441 1421830524E(—1) 

=  .73514100840602161576356511063802129E(— 2) 

=  .4571 9954438368161 04201 81 30444696542E(-5) 
=  .1227795 1 728998734846248338369272288E(— 1 ) 

=  .28782107488638388447298067656873723E(— 3) 

=  .73070314067824048001705945939654844E(— 6) 

=  .17832580880707298509927490339357089E(— 17) 

=  .44807240299498264499828954500884385E(— 4) 

=  .1000529345094521 1288323071593641690E(-8) 

=  .  188332416960281 64854032908626610601  E(-l 6) 
=  .7770581 09004299823785072765831 22027E(-52) 
=  .42842807196177999860791418123947654E(— 1) 

=  .11862011824909340596883700951819210 

=  .24732925056146632607609324152470353 

=  .34094562129787317503712966284227090 
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BRAT(A[2],  B[2],X[1]) 
BRAT(B[2],  A  [2],  1-X[2]) 
BRAT(B[2],  A  [2],  1-X[3]) 
BRAT(B[2],  A  [2],  1-X[4]) 
BRAT(B[3],  A[2J,  1-X[1]) 
BRAT(B[3],  A  [2],  1-X[2]) 
BRAT(B[3],  A[2],  1-X[3]) 
BRAT(B[3],  A  [2],  1-X[4]) 
BRAT(B[4],  A  [2],  1-X[1]) 
BRAT(B[4],  A  [2],  1-X[2]) 
BRAT(B[4],  A[2],  1-X[3]) 
BRAT(B[4],  A  [2],  1-X[4]) 
BRAT(A[3],  B[l],  X[l]) 
BRAT(A[3],  B[l],  X[2]) 
BRAT(A[3],  B[l],  X[3]) 
BRAT(A[3],  B[l],  X[4]) 
BRAT(A[3],  B[2],  X[l]) 
BRAT(A[3],  B[2],  X[2]) 
BRAT(A[3],  B[2],  X[3]) 
BRAT(B[2],  A  [3],  1-X[4]) 
BRAT(A[3],  B[3],X[1]) 
BRAT(B[3],  A[3],  1-X[2]) 
BRAT(B[3],  A[3],  1  -X  [3] ) 
BRAT(B[3],  A  [3],  1-X[4]) 
BRAT(B[4],  A  [3],  1-X[1]) 
BRAT(B[4],  A[3],  1-X[2]) 
BRAT(B[4],  A[3],  1  -X  [3] ) 
BRAT(B[4],  A[3],  1-X[4]) 
BRAT(A[4],  B[1],X[1]) 
BRAT(A[4],  B[l],  X[2]) 
BRAT(A[4],  B[l],  X[3]) 
BRAT(A[4],B[1],X[4]) 
BRAT(A[4],B[2],X[1]) 
BRAT(A[4],  B[2],  X[2]) 


=  .38001 1079761296945040771 20479083834 
=  .27117207094945604358130220571091526 
=  .60588505752236891712369301632731374E(— 1) 

=  .4540325883334931 1699534203586891 981  E(-4) 
=  .13405072148901327755227909052879592 
=  .5393515232200752080427 698 1 962034880E(  — 2) 
=  .19290796077732735634114110543000211E(— 4) 
=  .6008997 09201 53482181 8527955691 59478E(—1 6) 
=  . 1 264943809645527 066454072000999 1 680E  ( — 2) 
=  .52923753445854760109332516271631398E(— 7) 
=  .14588815980030198715137476274492143E(— 14) 

=  .78497514038461340403536969387757035E(— 50) 
=  ,32053821989521731672299996038770368E(— 5) 
=  .628200685374008430869961 84454438728E(-3) 

=  .17690074829443825633183038465251 117E(-1) 
=  .441 91 1 60942382014481 05704851 2463266 
=  .39791669199364443677769395264851494E(— 3) 
=  .40112235501938455812919777329719793E(— 1) 
=  .41006899961171290768594879542371616 
=  . 1 226460341 3482860075721 589941 537032E(— 2) 

=  ,35094060735746471970058683739185098E(— 1) 

=  .39702207630401328146312869741167449 
=  .95100701659741 858342243382133621 038E(— 2) 
=  . 1 25395621 1 352342546979905024 1820693E(— 12) 
=  .33218601813701179608821314194707665 
=  .48873314409464396177573572535759741  E(-3) 
=  .12965691544661661958807016151086044E(— 9) 
=  .34207572438829445357106628963957989E(— 44) 

=  .22686183750632833038237745473990532E(— 13) 
=  .37515536937999003719100338783203475E(— 7) 
=  .21915408366712445080427980539350023E(— 3) 

=  .29032021519120007197335268397917479 

=  .185651980281881 61 03373707081 1254010E(-10) 
=  .144898947665883394444240968391 48527E(— 4) 
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BRAT(A[4],B[2],X[3]) 
BRAT(B[2],  A  [4],  1-X[4]) 
BRAT(A[4],  B[3],  X[l]) 
BRAT(A[4],  B[3],  X[2]) 
BRAT(B[3],  A[4],  1-X[3]) 
BRAT(B[3],  A[4],  1-X[4]) 
BRAT(A[4],  B[4],  X[l]) 
BRAT(B[4],  A[4],  1-X[2]) 
BRAT(B[4],  A  [4],  1-X[3]) 
BRAT(B[4],  A  [4],  1-X[4]) 
BRAT(40.25, 45.25,  .25) 
BRAT(45.25, 40.25,  .50) 
BRAT(45.25, 40.25,  .25) 
BRAT(45.25, 40.25,. 15) 
BRAT(100.25, 201.25,  .325) 
BRAT(100.25, 201.25,  .33) 
BRAT(201.25, 100.25,  .665) 
BRAT(201.25, 100.25,  .66) 


=  .25011190127138313243352465595362702E(— 1) 

=  .10108292471982372254346080172398906E(— 1) 

=  .11 87258367812976464889702951 1 175989E(-6) 

=  .101798076475303081 12604293623106432E(—1) 

=  .33872352758566424339682566036491872 

=  .86292559297135160668356606705617935E(— 10) 

=  .23422666139670676884195587790023536E(— 2) 

=  .22696748842771075618260174519028746 
=  .64552803333196107540883097672076182E(— 5) 

=  .522489392890485 1 1 271 4021 560791 56800E(-38) 
=  .5737967625238081 1202923922414822492E(-5) 
=  .293378578061 85349044329069800003140 
=  .18931 66066778780701 926336541 8922239E( — 7 ) 
=  .  19972678050801 983781 084463824470771  E(- 1 5) 

=  .39598657181697412939230979816035992 

=  .4686459325768533402321 9947633433660 

=  .45800835607995502969459981579906802 

=  .38641541908994611 668339046999369105 


19.  (185— 256)— Page  89 

Call  DBRAT(a,  b,x,  y,  w,  wl,IERR)—  DOUBLE  PRECISION  BRATIO. 
Uses  same  arguments  as  BRATIO. 


20.  (257— 274)— Page  92 

Call  BSSLJ(z,n,w)-J  BESSEL  FUNCTION  COMPLEX  ARGUMENT,  INTEGER 
ORDER.. 

z=  (0, 1),  (5, 3),  (10,  -17.1),  n  =  0,  5,  10. 
z  =(0, 1) 

w(0)  =(1.2660658777520083355982446252147175,  0) 
w(5)  =(0,  .27146315595697187518107390515377734E(— 3)) 
w(10)=(-.27529480398368736252357102010027635E(-9),  0) 
z  =(5, 3) 

w(0)  =(-.82671066540395971142023564909882222,  3.2432199411359523046181066510440041) 
w(5)  =(.2233091 1878482278054181870553277152,  .98508733694458286449225062826498287) 
w(10)=(-.88641555041832623396338227470429441E(— 4),  — .86353649183906788526302676193032605E(— 2)) 
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z  =(10, 17.1) 

w(0)  =(-2293359.1837843503944023401281591159,  -727820.96263506654809560905605979146) 
w(5)  =(-823249.19677847131208538619627261896,  1112946.0824808045050358986659018367) 
w(10)=(-2898.958623354511941 7282522764953009,  263560.38741949811021174069631545538) 

21.  (275-292) 

Call  DBSLJ(z,n,  w)— J  BESSEL  FUNCTION  COMPLEX  ARGUMENT,  INTEGER 
ORDER.. 

Double  precision  version  of  BSSLJ  obtained  from  FCH.  Uses  same  arguments. 

22.  (293— 310)— Page  92 

Call  CBSSLJ(z,  u,  w)-J  BESSEL  FUNCTION  COMPLEX  ARGUMENT  AND 
ORDER.. 

=  (1,  1),  (5,3),  (10,  -17.25);  Zj  =  (0,1),  (5,3),  (10,-17.1) 
wi,j  =  w(i/j,  Zj);  i,j  =  l,...,3. 

wM  =  (.14642708861596860306081675249533466,  .52737913663990901653557434964427650E(-1)) 
w2>1  =  (.46218361676493820924597648673675168E(— 5),  ,25456792856760833967909202786165000E(-5)) 
w3,i  =  (-5575074.6647956050067554961459728768,  -5357286.9085164482330012058102014785) 
wi,2  =  (-.83227019009068576300571777487790209,  -.11619615843724546691699138790274513) 
w2>2  =  (.24429431173142058417192829176498112,  -,44252772655967690910905720945575648E(-1)) 
w3l2  =(7083286.9137785574760158909682753024,  4045803.5624379514501179491119759740) 
wi>3  =  (-2931371.2484970823912168697786184327,  10881912.336502746795574579146898874) 
w2>3  =(17579849.767255547205734221645520626, 125607208.16444575961851583496229498) 
w3j3  =  (.14990141503715764888724742102752180,  .62701690587612355833398868063692519E(-1)) 

23.  (31 1-328) -Page  91 

Call  DBSSLJ  (Z,  NU,  W)  -DOUBLE  PRECISION  VERSION  OF  CBSSLJ  . 

Z  =  (0,  1),  (5,  3),  (10,  -17.1),  NU  =  (1,  1),  (5,  3),  (10,  -17.25). 

24.  (329-346) -Page  99 

Call  BSSLK(MO,z,n,w)-K  BESSEL  FUNCTION  COMPLEX  ARGUMENT, 
INTEGER  ORDER. 

Same  arguments  for  n  as  BSSLJ,  z  =  (0,  4),  (5,  3),  (10,  —17.1),  MO  =  0. 

z  =(0,4) 

w(0)  =  (.26610451105001945410223079731868354E(— 1),  .62384146252142305380166549056420011) 
w(5)  = (-.20748123413740288540578776799012713,  -1.2501204889606838069891853734751489) 
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w(10)  =  (-280.12098217552229101545025349837889,  ,30636898683602146453207605247735934E(-3)) 
z  =(5, 3) 

w(0)  =  (— ,34121248887629061450927984359655133E(— 2),  .40776998107085526974196201704735308E(-3)) 
w(5)  =  ( — .  10382903044547709062774499608837364E (— 1 ) ,  ,16466569738005823397449971863723252E(-1)) 
w(10)  =  (2.5529708335454787244882116901948488,  .23322059884261042141739974174297953) 
z  =(10,  -17.1) 

w(0)  =  (,42104750178273550442610856630426844E(— 5),  -.12027131043135942524217995028976991E(-4)) 
w(5)  =  (.134962319359377350628301859671 88919E(— 4),  -,11451721272657599924162664202357980E(-4)) 
w(10) =  (.31396305911482406685020468010650756E(-4),  ,37919833264167302402012353819984313E(-4)) 

25.  (347-364) 

Call  DBSSLK(MO,z,n,w)-K  BESSEL  FUNCTION  COMPLEX  ARGUMENT, 
INTEGER  ORDER. 

Double  precision  version  of  BSSLK  obtained  from  FCH.  Uses  same  arguments. 

26.  (365— 382)— Page  99 

Call  CBESK(z,  u ,  w)— K  BESSEL  FUNCTION  COMPLEX  ARGUMENT  AND 
COMPLEX  ORDER. 

Uses  same  arguments  as  CBSSLJ—  wi;j  =  ;,  Zj). 

ww  =  (.40647504293828635018310643245122930,  -2.5249418710952837862596837388948793) 
w2)i  =  (8075.4464054492489165942412865748102,  -15508.059687733317145710512692794138) 
w3,i  =  (— .31358215995107745014533341187714286E(— 8),  -.97567321 1292285377921 67706290504644E(-9)) 
wi,2  =  (— .37126892365307062099384875555422001E(— 2),  -,77069030224451568902197301642716214E(-4)) 
w2,2  =  (— .88484788284610544522205958222540880E(-2),  -,29005444533438071634841325572526054E(-1)) 
w3,2  =  (.47411031953562064823534039929015534E(— 8),  -.94186612231164769089342007910597252E(-9)) 
w1>3  =  (.43341207234568979066457083646363663E(— 5),  -,11417790579579744437339079985538951E(-4)) 
w2,3  =  (.7271 66079901 76812148274009930 771 236E(— 5),  -,39571763205722002654735170788575282E(-5)) 
w3)3  =  (-. 114314831045567800687893733971 93914E(— 2),  -,86272020002795326178608532104303599E(-4)) 

27.  (383— 400)— Page  98 

Call  DBESK(Z,  NU,  W)-DOUBLE  PRECISION  VERSION  OF  CBESK. 

Uses  same  arguments  as  CBESK. 

28.  (401— 404)— Page  151 

x  =  ZEROIN(F, a, b,  AERR,RERR)— ZEROES  OF  CONTINUOUS  FUNCTIONS. 

F(u ,  r)  =  I (v ,  r)  -  K(u ,  r)  I,  K  BESSEL  FCNS  WITH  v ,  r  REAL.  Find  v  =  x  such 
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that  F(x,  r)  =  0  at  r  =  1.25,  3.5,  5.,  10.5  with  AERR  =  RERR  =  lOSPMPAR(l), 
a  =  0,  b  =  20. 

F(x,  1.25)  =0,  v  =  x  =  1.3569577601448969392878862075402274 
F(x,  3.5)  =0,  v  =  x  =  4.7859556026108308380903198411586656 

F(x,  5.0)  =0,  v  =  x  =  7.0547662294145426381757989398728754 

F(x,  10.5)  =0,  v  =x  =  15.360223427680893789116904053096803 

29.  (405— 408)— Page  151 

X  =  DZERO(F,  a,  b,  AERR,  RERR) -DOUBLE  PRECISION  VERSION  OF  ZEROIN. 
Same  input  arguments  as  ZEROIN  except  AERR=RERR=100DPMPAR(1),  Double 
precision:  DBSSLI  and  DBSSLK. 

30.  (409— 502)— Page  157 

Call  DRPOLY (A,  n,  ZR,  ZI,  NUM,  WK,  IWK)-DOUBLE  PRECISION  ROOTS  OF 
POLYNOMIALS. 

ZR(J&(J+1))  =»  ZR(J  +  1)  =  ZR(J),  ZI(J  +  1)  =  -ZI(J).  ZI  =  0,  if  not  listed. 

PI  =  3x4  -  16x3  -  3x2  +  13x  +  16 

A(l)  =  3,  A (2)  =  -16,  A(3)  =  -3,  A(4)  =  13,  A(5)  =  16. 

ZR(1&2)  =  -.66235897862237301298045442723904867, 

ZI(1&2)  =  ±.56227951206230124389918214490937306, 

ZR(3)  =  1.3247179572447460259609088544780973,  ZR(4)  =  16/3. 

P2  =  x3  -  lllx2  +  1130x  -  3000 

A(l)  =  1,  A(2)  =  -111,  A (3)  =  1130,  A (4)  =  -3000 

ZR(1)  =  5,  ZR(2)  =  6,  ZR(3)  =  100. 

P3  =  x2O+103x19+105x18+107x17+108x16+101°x15+10nx14+1012x13+1013x12+1014x11+ 
1016xlo+  1017x9+  1017x8+ 1018x7+  1019x6+  1019x5+  1019x4+  102Ox3+  102°x2  +  1019x+ 1019 


ZR(1) 

ZR(2&3) 

ZI(2&3) 

ZR(4) 

ZR(5) 

ZR(6&7) 

ZI(6&7) 

ZR(8) 


=  -901.23101299521587593226618428291990 
=  -49.8178965341880883285324369040581 18 
=  ±87.235177234039235432808002205227538 
=  -13.06444407624895220253180940339901 1 
=  -10.212353245230366224531415001988674 
=  -5.6146440952376396714047197455366775 
=  ±16.345595324801959495379204963055960 
=  -4.3896715520749770738308164416109363 
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ZR(9) 

ZR(10) 

ZR(11&12) 

ZI(11&12) 

ZR(13&14) 

ZI(13&14) 

ZR(15&16) 

ZI(15&16) 

ZR(17&18) 

ZI(17&18) 

ZR(19&20) 

ZI(19&20) 


=  -2.271 24346243023530692882271 77846490 

=  -1.1133698400337566260605048816290935 

=  -.4276593146381 1996386647708443938082E(-4) 

=  ±.31768590601815950904846274988714440 

=  +1.1 266346074981 447400220207122339536 
=  ±1 .926296567236464870248071 5208645235 

=  +2.25767002556886029431 67522772773882 

=  ±3.8726756637371263811836888992538821 

=  +5.86330080086494060001 1 6272590630581 
=  ±31.270717298930054322838041312035918 
=  +12.326025547042327860657919413394970 
=  ±9.5493492083624067302039056516068169 


P4  =  IIJ=1(x  -  k)  =  E“=o  A(j  +  l)xn-j,  n  =  20. 

A(l)  =  1,  A(2)  =  -210,  A (3)  =  20615,  A(4)  =  -1256850,  A(5)  =  53327946, 

A(6)  =  -1672280820,  A(7)  =  40171771630,  A(8)  =  -756111184500,  A(9)  =  11310276995381, 
A(10)  =  -135585182899530,  A(ll)  =  1307535010540395,  A(12)  =  -10142299865511450, 
A(13)  =  63030812099294896,  A(14)  =  -311333643161390640,  A(15)  =  1206647803780373360, 
A(16)  =  -3599979517947607200,  A(17)  =  8037811822645051776, 

A(18)  =  -12870931245150988800,  A(19)  =  13803759753640704000, 

A(20)  =  -8752948036761600000,  A(21)  =  2432902008176640000 
ZR(k)  =  k,  ZI(k)  =  0,  k  =  1, . . .  ,  20. 


31.  (503— 514)— Page  157 

Call  DCPOLY(AR,  AI,  n,  ZR,  ZI,  NUM,  DWK) -DRPOLY,  EXCEPT  COEFFICIENTS 
ARE  COMPLEX. 

PC1=  (Z  +  A)6,  A=1+7.3I. 

AR(1)  =  1,  AR(2)  =  6,  AI(2)  =  43.8,  AR(3)  =  -784.35,  AI(3)  =  219, 

AR(4)  =  -3177.4,  AI(4)  =  -7342.34,  AR(5)  =  37816.2615,  AI(5)  =  -22903.020, 

AR(6)  =  82003.3230,  AI(6)  =  101262.27558,  AR(7)  =  -109535.214789,  AI(7)  =  116647.75558 
ZR(k)  =  -1,  ZI(k)  =  -7.3,  k  =  1, ...  ,6. 

32.  (51 5-524) -Page  175 

Call  S AXPY (n,  a,  X,  kx,  Y,  ky) -VECTOR  ADDITION. 

n  =  5,  a  =100,  kx  =  2,  ky  = -2,  X  =  (1,0, 2, 0,3, 0,4, 0,5,0),  Y=  (11,0, 10,9,0,8,0,7,0) 
Output  :Y  =  (511,0,410,0,309,0,208,0,107,0) 
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33.  (525-534) 

Call  DAXP Y (n,  a,  X,  kx,  Y,  ky)-VECTOR  ADDITION. 

Double  precision  version  of  SAXPY  obtained  from  FCH.  Uses  same  arguments. 

34.  (535-634) -Page  216 

Call  KROUT(MO,  n,  m,  A,  ka,  B,  kb,  IERR,  INDEX,  TEMP)-INVERT  MATRIX 
MO  =  0,  n  =  10,  m  =  — 1,  A  =  Hilbert  matrix  of  order  10,  ka  =  10. 

Inverse  stored  in  A.  Elements  of  inverse  given  in  APPENDIX  A. 

35.  (635-734) 

Call  DKROUT(MO,  n,  m,  A,  ka,  B,  kb,  IERR,  INDEX,  TEMP) -INVERT  MATRIX. 
Double  precision  version  of  KROUT  obtained  from  FCH.  Uses  same  arguments. 

36.  (735-835) -Page  217 

Call  MSLV (MO,  n,  m,  A,  ka,  B,  kb,  D,  RCOND,  IERR,  IWK,  WK). 

Same  arguments  as  KROUT.  Output:  Inverse  of  A  stored  in  A. 

D  =  Determinant  of  A  =  .21641792264314918690605949836507259E(-52) 

37.  (836-936) -Page  21 7 

Call  DMSLV(MO,  n,  m,  A,  ka,  B,  kb,  D,  RCOND,  IERR,  IWK,  WK) 

Double  precision  version  of  MSLV.  Same  arguments  as  KROUT.  Output:  Inverse  of  A 
stored  in  A. 


38.  (937-942) -Page  221 

Call  SLVMP(MO,  n,  A,  ka,  b,  X,  WK,  IWK,  IND)-SOLUTION  OF  REAL  EQUATIONS 
WITH  ITERATIVE  IMPROVEMENT. 

MO  =  0,  n  =  3,  ka  =  3.  Call  SLVMP  again  with  MO  =  1. 

SLVMP  solves  AX  =  b  =  bj  with  MO  =  0.  Then  another  call  is  made  to  SLVMP  with 
MO=l  and  it  solves  AX  =  b  =  b2. 


A= 


1  5  2  \ 
6  3  7 
4  8  5  / 


/ 


7 T 


exp(-l) 

ln(2) 


X  (with  bi)  =  4/3,  1/3,  -1 

X  (with  b2)  =  -5.3340450253088614142013505296964594, 

-.18671846323305651879169987951202407,  4.7046149975319686233112466552680414 
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39.  (943-948) 

Call  DSLVMP(MO,  n,  A,  ka,  b,  X,  WK,  IWK,  IND) 

SOLUTION  OF  REAL  EQUATIONS  WITH  ITERATIVE  IMPROVEMENT. 

Double  precision  version  of  SLVMP  obtained  from  FCH.  Uses  same  arguments. 

40.  (949-983) -Page  353 

Call  EIGV(IBAL,  A,  ka,  n,  WR,  WI,  ZR,  ZI,  IERR) 

FINDS  EIGENVALUES  AND  EIGENVECTORS  OF  A  MATRIX. 

IBAL=1,  Input  matrix  A,  ka=No.  of  rows  of  A,  n=order  of  A,  (ZR(K,J),  ZI(K,J))=the 
Kth  eigenvector  of  the  Jth  eigenvalue  of  A.  Output  requires  IERR=0.  Let  H(5)  denote 
the  Hilbert  matrix  of  order  5  (n=5).  First  the  eigenvalues  of  A  =  H(5)  x  H(5)  are  given. 
Then  the  eigenvalues  and  eigenvectors  for  A=H(5)  are  listed. 

A  =  H(5)  x  H(5) 

WR(1)  =  . 1 08104756108755743071 08389406252299E(— 1 0) 

WR(2)  =  .93573610968340105856166349581 076979E(-7) 

WR(3)  =  .13013086513839305374962217453885253E(— 3) 

WR(4)  =  .43486520331705900754379338104236462E(— 1) 

WR(5)  =  2.4556478684714427530446185082147609 

A  =  H(5) 

WR(1)  =  .32879287721 718629571 150047605447314E(-5) 

ZR(1, 1)  =  .12197366118335218941776511729954141E(— 1) 

ZR(2, 1)  =  -.23054351082321566090721878767252246,  ZR(3, 1)  =  1 
ZR(4, 1)  =  -1.5156978981503881831443663648351046 
ZR(5, 1)  =  .74332785300803664576594824839796720 
WR(2)  =  .305898040151 19172687949784069272283E(-3) 

ZR(1,2)  =  .7067 02262 1 08752454469500 1 1 9344285 1 9E( — 1 ) 

ZR(2, 2)  =  -.64833602593662613024121970477768598,  ZR(3, 2)  =  1 
ZR(4, 2)  =  .34917863233062409987909081721800737 
ZR(5, 2)  =  -.83554293387428295255819834635554228 
WR(3)  =  .11407491623419806559451458866589345E(— 1) 

ZR(1,3)  =  -1.7783959618547884961815728224945915 

ZR(2, 3)  =  6.0114771493213747655493913824738878,  ZR(3, 3)  =  1 

ZR(4, 3)  =  -2.5700750831853313392911158332983055 

ZR(5, 3)  =  -4.6922210619044074735319696492611082 


16 


NSWCDD/TR-98/75 


WR(4)  =  .20853421861101333590500251006882006 
ZR(1 , 4)  =  -1.4165794178621925206014926444904337 
ZR(2, 4)  =  .64939656106525325387143805590856927,  ZR(3, 4)  =  1 
ZR(4, 4)  =  1.0447810384175732744990567015288115 
ZR(5, 4)  =  1.0097363119375138848812509913358092 
WR(5)  =  1.5670506910982307955330110055207246 
ZR(1,5) =  2.3877691630491462137700003189820182 
ZR(2, 5)  =  1.3862597946271820256958630223217197,  ZR(3,5)  =  1 
ZR(4, 5)  =  .78810960688381844435363097568288750 
ZR(5, 5)  =  .65247760861091767322074653366065188 

41.  (984- 101 8) -Page  357 

Call  DEIGV(IBAL,  A,  ka,  n,  WR,  WI,  ZR,  ZI,  IERR) -FINDS  DOUBLE  PRECISION 
EIGENVALUES  AND  EIGENVECTORS  OF  A  MATRIX. 

Same  input  as  EIGV. 


42.  (1019—1 026 )  —Page  377 

Call  LLSQ(m,  n,  A,  ka,  B,  kb,  L,  WK,  IWK,  IERR)— LEAST  SQUARES  (LSQ) 
SOLUTION  OF  A  LINEAR  SYSTEM  OF  EQUATIONS. 

Rank  of  input  matrix  A=n.  L=No.  of  LSQ  problems  (No.  of  columns  of  B).  WK  and 
IWK=work  space.  Output  X  stored  in  B.  Norm  stored  in  B(n+1,  J),  J=  1, . . .  ,  L. 
Output  IERR  must  be  0. 


A  =  Ai  = 


A  =  A2  = 


/l 
1 
1 

V1 

/ 1 
1 
1 

v1 


1  0 

(  3  \ 

2  4 

B  = 

10 

0  0 

3 

-1  1  j 

CO 

►-* 

o 

1 ) 

3  \ 

2 

B  = 

10 

0 

3 

-ij 

9/10 ; 

X  =  Xi 


327/200  \ 
301/200 
V  49/40  / 


X  =  X2  = 


143/50  \ 
273/100  ) 
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a  =  a3  = 

(  1  2  3  ^ 
4  5  6 

B  = 

3  1 
10 

x  =  x3  = 

/  -119/15  \ 
499/60 

7  8  10 

^  10  11  12  ) 

3 

9/10  ) 

-4/5  ) 

||  (A3X3-B)  ||norm=  46.705 

43.  (1027— 1034)— Page  385 

Call  DLLSQ(m,  n,  A,  ka,  B,  kb,  L,  WK,  IWK,  IERR)  —DOUBLE  PRECISION  LEAST 
SQUARES  (LSQ)  SOLUTION  OF  A  LINEAR  SYSTEM  OF  EQUATIONS. 

Same  input  as  LLSQ. 

44.  (1035-1037) -Page  409 

Call  OPTF(F,  n,  RERR,  ITER,  X,  FVAL,  IND,  WK) 

FINDS  LOCAL  MINIMUM  OF  F  AT  X=  (X(l),  X(2),  . . .  ,  X(n)). 

RERR=  Relative  accuracy  of  F,  FVAL=F(X),  ITER=No.  of  iterations,  WK=work  space. 
Output:  IND=1  OR  2.  F  obtained  from  user  subroutine  Call  F(n,X,  FVAL). 

F  =  Fi  =  x2  —  2x  +  y2  +  2y,  Initial  guess  =  X(l)  =  —5,  X(2)  =  5. 

Output :  X(l)  =  1,  X(2)  =  -1,  FVAL  =  -2. 

F  =  F2  =  x3  +  y2  —  3x,  Initial  guess  =  X(l)  =  X(2)  =  0. 

Output :  X(l)  =  1,  X(2)  =  0,  FVAL  =  -2. 

F  =  F3  =  2x2  +  2y2  +  z2  —  2xy  —  2yz  —  4x  +  4y  —  2z  +  3, 

Initial  guess  =  X(l)  =  X(2)  =  X(3)  =  0. 

Output :  X(l)  =  1,  X(2)  =  0,  X(3)  =  1,  FVAL  =  0. 

45.  (1038-1040) 

Call  DOPTF(F,  n,  RERR,  ITER,  X,  FVAL,  IND,  WK) 

Double  precision  version  of  OPTF  obtained  from  FCH.  Uses  same  arguments. 

46.  (1041 -1059)-Page 413 

Call  SMPLX(A,  B,  C,  ka,  m,  n,  IND,  IB,  X,  z,  ITER,  MIT,  NLE,  NGE,  BI,  WK,  IWK) 
Find  X  when  z=C*X  is  maximized  subject  to  A*X(<,  >,  =)B.  DIM(A)  =  m  x  n, 
DIM(B)  =  m,  DIM(C)=n.  IND  =  0  SMPLX  selects  starting  input,  IB  not  used, 
ITER=No.  of  iterations,  MIT=Max  No.  of  ITER’s  allowed.  NLE(NGE)=  No.  of  <  (>). 
Output:  X={xj},  z.  Four  CASES 
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A  =  A-i  — 


-X(l)  +  X(2)  <  5 
X(l)  +  4X(2)  <  45 

2X(1)  +  X(2)  <  27 

3X(1)  —  4X(2)  <  24 


(-1 

1  \ 

(  <h  ^ 

1 

4 

<45 

B  —  Bx  — 

2 

1 

<27 

3 

~4J 

U24^ 

C  =  C!  = 


NLE  =  4,  z  =  MAX{X(1)  +  X(2)}  =  9  +  9  =  18. 


A  =  Ao  = 


(  5  2  4  \ 
1  1  2 
4  3  6 
1  1  1 


B  =  B2  = 


(  <20  \ 
<  8 
>5 
=  6 


C  =  C2 


2 

1 

V-3 


z=MAX{2X(l)  +  X(2)  -  3X(3)}  =  2x8/3+  10/3  +  0  =  26/3. 

A  =  A3  = 

z  =  MAX{2X(1)  +  ,1X(2)  +  10X(3)}  =  2  x  4/3  +  .1  x  80/3  +  0  =  16/3. 


200  -10 

2000  \ 

(  <0  \ 

(  2 

10 

1 

100 

B  =  B3  =  <40 

c  =  c3  = 

.1 

-1 

.2 

10  j 

\  >  !  ) 

10 

A  =  A4  = 


/  0  0  3  0  1  1  0  \ 

0  12-10  0  0 
1  0  0  0  0  -1  0 

0  0  1  0  0  1  1 


B  =  B4  = 


/  =  2  \ 
=  20 
=  10 
=  3 


C  =  C4  = 


z  =  MAX{X(1)  +  X(2)  +  X(3)  -  3X(4)  +  X(5)  -  X(6)  +  3X(7)}  = 

10+  20+  0-  3x0+  2-  0+  3x3=  41. 


/  1  \ 
1 
1 

-3 
1 

-1 

\  3  ) 


47.  (1060-1078) 

Call  DSMPLX(A,  B,  C,  ka,  m,  n,  IND,  IB,  X,  z,  ITER,  MIT,  NLE,  NGE,  BI,  WK,  IWK) 
Double  precision  version  of  SMPLX  using  FCH.  Uses  same  arguments. 

48.  (1079—1086)  —Page  413 

Call  SSPLX(TA,  ITA,  JTA,  B,  C,  m,  n,  IND,  IB,  X,  z,  ITER,  MIT,  NLE,  NGE,  BI,  WK,  IWK) 
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Generates  the  same  results  as  SMPLX.  Used  when  A  is  sparse.  TA,  ITA,  JTA  specify 
AT  (TRANSPOSE)  in  sparse  format— Page295.  For  CASE  4  of  SMPLX  above 
TA  =  TA4=(1,  1,  1,  1,  1,  3,  2,  1,  1,  -1,  1,  1,  1,  1,  -1,  1,  1,  1)T 
JTA  =  JTA4  =  (1,  10,  2,  9,  3,  8,  9,  11,  4,  9,  5,  8,  6,  8,  10,  11,  7,  11)T 
ITA  =  ITA4  =  (1,  3,  5,  9,  11,  13,  17,  19)T,  m  =  11,  n  =  7 
B  =  B4  =  (0,  0,  0,  0,  0,  0,  0,  2,  20,  10,  3)T,  NLE  =  0,  NGE  =  7 

49.  (1087-1094) 

Call  DSSPLX(TA,  ITA,  JTA,  B,  C,  m,  n,  IND,  IB,  X,  z,  ITER,  MIT,  NLE,  NGE,  BI,  WK,  IWK) 
Double  precision  version  of  SSPLX  obtained  by  FCH.  Uses  same  arguments. 

50.  (1095— 1 099)  —Page  433 

Call  CHEB Y (a,  b,  F,  G,  PHI,  EPS,  ITER,  MIT,  L,  m,  P,  Q,  ERR,  IERR,  WK) 

MINIMAXFIT. 

CHEBY  is  a  double  precision  routine,  a  <  X  <  b,  EPS  =  10-4,  MIT  =  40,  WK=work 
space,  G=F(X),  PHI= Argument  of  polynomials.  Output:  P,  Q;  ERR  specifies  relative 
agreement  between  F  AND  P/Q  at  critical  points.  IERR  should  be  0. 

CASE:  Fit  f(X,K),  THE  ELLIPTIC  INTEGRAL  OF  THE  FIRST  KIND  (Page  107), 
where  K=.9.  PHI=  X2,  F(X)=f(X,K)/X  on  [a  =  0,  b  =  tt/2].  The  final  result  X  x  P/Q 
is  compared  with  f(X,K)  at  X  =  7r/10,  tt/5,  7t/4,  27t/5,  tt/2.  L=5  and  m=6.  The  elliptic 
integral  is  evaluated  by  Call  DELLPI  —Page  109. 

f(Xj,  K)  =  .31841393498182425425662894958775021 
f(X2,  K)  =  .66408837881648915721400693335979825 
f(X3,  K)  =  .85794019788551098164818816682426150 
f(Xi,K)  =  1.6019861670355365511633360018405787 
f(X5,  K)  =  2.2805491384227702046137519445555304 

51.  (1100—1104) —Page  447 

Call  HTRP(n,  X,  Y,  A,  WK,  IERR) -HERMITE  INTERPOLATION. 
n=10  for  case  below.  Xj  and  the  associated  function  and  the  first  derivative  values  yjt  yj 
are  stored  in  X  =  (X(2j  —  1)  =  Xj,  X(2j)  =  xj ),  Y  =  (Y(2j  —  1)  =  yj,  Y(2j)  =  yj ), 
j  =  l,  •••,5. 

CASE  :  y  =  EXP(x),  Xj  =  .1,  .2,  .3,  .4,  .5. 


20 


NSWCDD  /TR-98 /75 


yi  =  yi  =  EXP(.l)  =  1.1051709180756476248117078264902467 
y2  =  y>2  =  EXP  (.2)  =  1.2214027581601698339210719946396742 
y3  =  /3  =  EXP(.3)  =  1.3498588075760031039837443133280073 
y4  =  /4  =  EXP  (.4)  =  1.4918246976412703178248529528372223 
y5  =  =  EXP  (.5)  =  1.6487212707001281468486507878141636 

The  interpolation  polynomial  obtained  from  HTRP  with  its  coefficients  stored  in  A  is 
converted  by  using  the  routine  PCOEFF(c>!,n,X,A,C,T)  —Page  449,  (a  =  0,  T=work 
space)  from  A(l)  +  A0  +  l)(x  -  xx)  •  •  •  (x  -  xj)  to  C0  +  !)xj-  The  latter  is 
used  for  five  arguments  and  compared  with  the  MAPLE  values  given  below: 

EXP(.15)  =  1.1618342427282831226166202143316552 
EXP(.25)  =  1.2840254166877414840734205680624365 
EXP(.35)  =  1.4190675485932572482703956619398724 
EXP(.45)  =  1.5683121854901688111795997746932154 
EXP(.55)  =  1.7332530178673952368219167671373288 

52.  (1105-1109) 

Call  DHTRP(n,  X,  Y,  A,  WK,  IERR)  -HERMITE  INTERPOLATION. 

Double  precision  version  of  HTRP  obtained  from  FCH.  Uses  same  arguments. 

53.  (1110— 1127)— Page 451 

Call  PFIT(n,  m,  X,  Y,  A,  RNORM,  PHI,  WK,  IERR)-LSQ  POLYNOMIAL  FIT. 

Given  (x*,  y;),  i  =  1,  ...  ,  m,  PFIT  fits  a  unique  polynomial  of  degree  n  by  minimizing 
RNORM2  =  £r=i[P(xO  —  yj]2.  The  polynomial  coefficients  are  stored  in  A  with 
P(X)  =  E“=0A(j  +  1)XL  PHI  and  WK  are  work  spaces. 

CASE  :  n  =  6,  m  =  10,  y  =  LOG(x),  X(i)  =  i/2,  Y(i)  =  LOG(X(i)),  i  =  1,  . . .  ,  m 
X(l)  =  .5,  Y(l)  =  -.69314718055994530941723212145817657 

X(2)  =  1.0,  Y(2)  =  0.0 

X(3)  =  1.5,  Y(3)  =  .40546510810816438197801311546434914 

X(4)  =  2.0,  Y(4)  =  .69314718055994530941723212145817657 

X(5)  =  2.5,  Y(5)  =  .91629073187415506518352721176801107 

X(6)  =  3.0,  Y(6)  =  1.0986122886681096913952452369225257 

X(7)  =  3.5,  Y(7)  =  1.2527629684953679956881206219850032 

X(8)  =  4.0,  Y(8)  =  1.3862943611198906188344642429163531 

X(9)  =  4.5,  Y(9)  =  1.5040773967762740733732583523868748 

X(10)  =  5.0,  Y(10)  =  1.6094379124341003746007593332261876 
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A(l)  =  -1.9800785652725251593569509944240339 
A(2)  =  3.4868393335634220135011658377975746 
A(3)  =  -2.2320116785443095629139842118019607 
A  (4)  =  .91658519020745719115335408583374835 
A(5)  =  -.21880135089709283048347902999887716 
A(6)  =  .27681 436863 1 4478666471 4633861 698448E(—1 ) 
A(7)  =  — .14317677628543903352379990489323698E(— 2) 
RNORM  =  .404910516290816676155263681 18077578E(-2) 


54.  (1128-1145) 

Call  DPFIT(n,  m,  X,  Y,  A,  RNORM,  PHI,  WK,  IERR) 

Double  precision  version  of  PFIT  obtained  from  FCH.  Uses  same  arguments. 


55.  (1146— 1179)— Page 455 

Call  CBSPL(X,  Y,  A,  B,  C,  n,  I,  J,  W,  V,  IERR)-CUBIC  SPLINE. 

Input  values:  X  =  (xi, . . .  ,  xn),  Y  =  (yi, . . .  ,  yn). 1  =  J  =  2,  W  =  V  =  0. 

Output:  Arrays  A,  B,  C  which  contain  the  coefficients  of  the  spline  F,  (x)  where 

Fj(x)  =  yj  +  a;(x  -  x5)  +  b;(x  -  Xj)2  +  Cj(x  -  Xi)3,  i  =  1, . . .  ,  n  —  1. 

The  Maple  result  gives  Fj(x)  =  y~]|_n  djj  xL  Hence  for  comparison  one  has 
dito  =  ((-Ci  x5  +  bi)  x5  -  ai)  xj  +  y} 
di,i  =  (3ciXj  -2bi)xs  +  ai 
di,  2  =  b;  3  Ci  X; 
d\,3  =  Cj 

Maple  results  for  fitting  y=f(x)  on  [0,2],  where  f(x)  =  sinh(x)/H  —  x3/12, 


H  =  sinh(2)  =  3.6268604078470187676682139828012617,  n  =  9,  are  : 


yi,  =  y(xi  =  0.00),  = 

Y2,  =  y(x2  =  .250),  = 

y3,  =  y(x3  =  .500),  = 

y4,  =  y(x4  =  -750),  = 

ys,  =  y(x5  =  1.00),  = 

ye,  =  y(x6  =  1-25),  = 

y7,  =  y(x?  =  1.50),  = 

ys  =  y(x8,=  1.75),  = 

y9,  =  y(x9  =  2.00),  = 


0 

.6834832732532345651 7151 982082545372E(-1 ) 

.13326002526399426077516293656568618 

.19157338375063409925452199590007833 

.24069380349860936645415534327974183 

.27892161687257251424544092568240252 

.30583613391569788177127850212405712 

.32276004168403583070008267835042497 

.33333333333333333333333333333333333 
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di,o  =  0 

dM  =  .27571464554407011163764157784986716 
di,2  =  0 

d1>3  =  — .37141379884420569104538392314970775E(— 1),  x  <  .25 

d2,0  =  — .45374977511776094409446680123941891E(— 4) 

d2,i  =  .27625914527421142477055493801135446 

d2)2  =  — .21779989205652525316534406459492108E(— 2) 

d2)3  =  — .34237381323666899062333804787038494E(— 1)),  x  <  .5 

d3)0  =  — .79273026732965099034843331664170571E(— 3) 

d3ii  =  .28074327701311867414618885783046105 

d3)2  =  —.11 146262398379751282921 2802841 62377E(— 1 ) 

d3)3  =  —.2825853900512389989482191 1 694896383E(— 1 ) ,  x  <  .75 

d4i0  =  —  .47 90300488077 50184391 742343 1 1 353492E  ( — 2) 

d4ii  =  .29673355789611007756046481828843562 

d4;2  =  — .32466636909034955835289227561461808E(— 1) 

d4>3  =  —.18782817 0003882534271 0282401 6096636E  ( — 1 ) ,  x  <  1. 

d5)0  =  —.181 84195832899588642442061 780969739E(— 1 ) 

d5,i  =  .33691524393057633795603873333793879 

d5>2  =  —  .7264832294350 121 623086314261 0964978E(— 1 ) 

d5>3  =  — .53889216555661666285781856662622460E(— 2),  x  <  1.25 

d6,0  =  —  .5490932872621 3634961 930986334089364E(— 1 ) 

d6jl  =  .42505556287453004912281215226542589 

d6)2  =  -.14316057809866418516428187775295466 

d6i3  =  .  1 341 43463858 1 0625087000143704935002E(— 1 ) ,  x  <  1.5 

d7,0  =  -.13243024816007326936780977901232400 

d7|1  =  .58009740174224931793456973762189515 

d7)2  =  -.24652180401047703103878693465726750 

d7,3  =  .36383507699546813059112378572560075E(— 1),  x  <  1.75 

d8)0  =  -.33409108003353744375636485161886800 

d8ll  =  .92580168495390218831494986209025630 

d8>2  =  -.44406710870285009982757557721061673 

d8)3  =  .7401 1 184783808349971262596201769454E(— 1 ),  x  >  1.75 
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A  comparison  is  also  made  between  the  generated  spline  F;  and  f(x)  at  x— .1  and  1.6. 
Maple  gives: 

Fj(.l)  =  .2753432317452259059465961 9392671745E(-1) 
f(.l)  =  .27534699553492088299860036105729450E(— 1) 

F7(1.6)  =  .31365662389804818615833155109330953 
f(l  .6)  =  .31365960437678305898304387504099094 

56.  (1180-1213) 

Call  DCBSPL(X,  Y,  A,  B,  C,  n,  I,  J,  W,  V,  IERR)— CUBIC  SPLINE. 

Double  precision  version  of  CBSPL  obtained  from  FCH.  Uses  same  arguments. 

57.  (1214— 1214)— Page  465 

IntF  =  CSINT(X,  Y,  A,  B,  C,  N,  a,  (3) -INTEGRAL  OF  A  SPLINE  . 

X,  Y,  A,  B,  C,  N,  f(x),  F(x)  are  given  above  in  CBSPL.  Considering  the  integral  of 
f(x)(F(x))  from  0  to  2  requires  a  =  0,  /3  -  2.  It  is  denoted  by  Intf(0,2)(IntF(0,2)). 
Thus  Intf(0,2)=  (COSH(2)  -  1)/H  -  1/3.  A  comparison  is  made  with 
lntf(0, 2)  =  .42826082262243155478612494927146026  True  value 

IntF(0, 2)  =  .42825482879644533223662485594509265  Maple  spline  result. 

58.  (1215-1215) 

DIntF  =  DCSINT(X,  Y,  A,  B,  C,  N,  a,  /?) -INTEGRAL  OF  A  SPLINE  . 

Double  precision  version  of  CSINT  obtained  from  FCH.  Uses  same  arguments. 

59.  (1216— 1221)— Page 525 

Call  SFVAL(n,  X,  Y,  Z,  m,  XI,  YI,  ZI,  IADJ,  IEND,  DZ,  IERR) 

SURFACE  INTERPOLATION. 

Given  surface  points  (xj,  y,),  i  =  1,  ...  ,  n,  stored  in  X  and  Y;  a  triangular  mesh  is 
generated  by  TRMESH(n,X,Y, IADJ, IEND, IERR) -Page  523  which  stores  the  results 
in  IADJ  and  IEND.  Let  f(x,y)  specify  a  surface.  The  numerical  partial  derivatives 
Fx(xj,  yj)  and  Fy(x;,  y,)  required  by  SFVAL,  are  generated  by  GRADL(n,X,Y,Z,IADJ, 
IEND, DZ, IERR) -Page  524.  Its  results  are  stored  in  DZ(k,n),  where  the  first  row  (k=l) 
contains  the  partials  with  respect  to  x  and  the  second  row  (k=2)  the  partials  with  respect 
to  y.  The  interpolation  is  to  be  carried  out  at  the  points  (xj,  yj),  j  =  1,  . . .  ,  m  stored 
in  XI  AND  YI.  The  interpolated  values  of  f,  F(xj,  yj)  are  stored  in  ZI. 

Let  f(x,y)=Circular  coverage  function,  evaluated  by  Call  CIRCV(R,  D,  I)  —  Pagel25 
with  I  =  1; Xj  =  R,  y,  =  D;  R,  D  =  (i  —  l)/8,  i  =  1, . . .  ,  n  =  80.  (6400  points) 
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XI  =  (.0625,  7.5625,  1.0625,  2.0625,  5.0625,  8.0625),  YI  =  (0,  9.0625,  1.0625,  2.0625, 
4.0625,  3.0625).  The  values  to  compare  are  given  by: 
f(.0625,  0)  =  .19512188925245272899573409173264047E(— 2) 

F(.0625,  0)  =  .1997269391 36652E(-2) 

f(7.5625, 9.0625)  =  .59324966894443808464458766563315521E(-1) 

F(7.5625, 9.0626)  =  .593298922404450E(-1) 

f(l  .0625,  1.0625)  =  .28253832617966374058497114901082367 

F(1.0625,  1.0625)  =  .282556449740250 

f(2.0625,  2.0625)  =  .39990590720072733397670102612204638 

F(2.0625,  2.0625)  =  .399917529883358 

f(5.0625, 4.0625)  =  .81301989279372316483739066020196562 

F(5.0625, 4.0625)  =  .812995939477743 

f(8.0625, 3.0625)  =  .99999952719147575100918544926929062 

F(8.0625,  3.0625)  =  .999999539276169 

60.  (1222-1227) 

Call  DSFVAL(n,  X,  Y,  Z,  m,  XI,  YI,  ZI,  IADJ,  IEND,  DZ,  IERR) 

SURFACE  INTERPOLATION. 

Double  precision  version  of  SFVAL  obtained  from  FCH.  Uses  same  arguments. 

61.  (1228— 1233)— Page 531 

Call  QAGS(f,  a,  b,  AERR,  RERR,  z,  ERROR,  NUM,  IERR,  L,  m,  n,  IWK,  WK) 
EVALUATION  OF  INTEGRALS. 

z  contains  the  output  value  for  the  estimate  of  the  integral  of  f  from  a  to  b.  ERROR=the 
estimated  absolute  error  in  z;  NUM  =  no.  of  points  at  which  f  was  evaluated.  IERR 
should  be  zero.  IWK(L)  and  WK(m)  are  work  spaces.  n=  no.  of  subintervals  used. 
CASE:  Circular  coverage  function  tested  with  the  input  taken  from  XI  and  YI  of  SFVAL 
above.  AERR=RERR=50  SPMPAR(l).  Uses  BESI-Page  97. 

62.  (1234— 1239)— Page 534 

Call  DQAGS(F,  a,  b,  AERR,  RERR,  z,  ERROR,  NUM,  IERR,  L,  m,  n,  IWK,  WK) 
EVALUATION  OF  INTEGRALS. 

Double  precision  version  of  QAGS  with  same  CASE.  AERR=RERR=50DPMPAR(1). 
Uses  double  precision  version  of  BESI  obtained  from  FCH. 

63.  (1240-1245) -Page  533 

x  =  QSUBA(F, a, b, RERR, MCOUNT, ERROR, IERR) 
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EVALUATION  OF  INTEGRALS  . 

Variables  that  also  appear  in  QAGS  have  the  same  meaning,  except  for  ERROR  which 
specifies  the  relative  error  here  and  RERR  =  max[50SPMPAR(l),  1E(  14)].  Input 
MCOUNT  specifies  the  maximum  no.  of  times  F  is  to  be  evaluated.  CASE:  Same 
as  used  above  in  QAGS. 

64.  (1246-1251) 

X  =  DQSUBA(F,  a,  b,  RERR,  MCOUNT,  ERROR,  IERR) 

EVALUATION  OF  INTEGRALS. 

Double  precision  version  of  QSUBA  obtained  from  FCH.  Uses  same  arguments.  Be¬ 
cause  of  the  limited  accuracy  of  the  modified  Bessel  function  Io  which  appears  in  F, 
RERR=max[50DPMPAR(l),  IE-14]. 

65.  (1252— 1252)— Page  539 

Call  QAGI(F,  a,  MO,  AERR,  RERR,  z,  ERROR,  NUM,  IERR,  L,  m,  n,  IWK,  WK) 
EVALUATION  OF  INTEGRALS  OVER  INFINITE  INTERVALS. 

For  MO=  -1,  z  =  F(x)dx;  for  MO  =  1,  z  =  /a°°  F(x)dx;  for  MO  =  2,  z  =  F(x)dx. 
Except  for  a,  variables  that  also  appear  in  QAGS  have  the  same  meaning. 

CASE:  With  b  =  \/8, 

/OO  /7F  "  1 

e"t2  Y2/3(bt)dt  =  tan(7r/3)I1/3(l)  +  7rC0S(7r/3)Kl/3(1)j  > 

see  [3,  Page  487,  Eq.  (11.4.30)].  Hence  a=0,  MO  =  1,  AERR=RERR=  1 00 SPMPAR(  1 ) . 

F  gives  the  integrand.  Uses  Bessel  function  Y  by  Call  CBSSLY— page  95. 

ZTRUE  =  -.69218599490499706343637364910052890  (From  MAPLE). 

66.  (1253— 1253)— Page 540 

Call  DQAGI(F,  a,  MO,  AERR,  RERR,  z,  ERROR,  NUM,  IERR,  L,  m,  n,  IWK,  WK) 
Double  Precision  version  of  QAGI;  same  input  is  used  with  DBSSLY. 

67.  (1254— 1263)— Page  551 

Call  ODE(F,  n,  Y,  T,  TOUT,  RERR,  AERR,  IND,  WK,  IWK) 

ADAPTIVE  ADAMS  DE  INITIAL  VALUE  SOLVER 

ODE  solves  the  set  of  n  DE’s  y'(t)  =  f(t,y(t))  numerically  to  obtain  y(t)  on  T=t=a  to 
TOUT=t=b,  b  >  a,  with  a  set  of  given  initial  conditions  y(a).  Y  contains  (yi(t), . . .  ,  yn(t))- 
F  computes  (yi(t), . . .  ,y/n(t))  using  the  DE’s  and  subroutine  F(t,Y,DY).  Two  CASES 
are  carried  out.  RERR  =  AERR  =  10SPMPAR(1),  IND  =  -1. 
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CASE1: y y"  —  (y/)2  — y2 log(y)  =  0;  y(0)  =  y'(0)  =  1;DE  =*  y"(0)  =  1,  [4,p.  571,6.113]. 
Hence  for  yi(0)  =  Y(l)  =  1,  y[(0)  =  y2(0)  =  Y(2)  =  1,  DY(2)  =  Y(2)2/Y(l)  + 
Y(l)  log(Y(l).  The  solution  of  the  DE’s  is  y(t)  =  esmh<^,  y'(t)  =  y(t)  cosh(t). 

Output(l) :  With  a  =  0,  b  =  1  : 

Y(l)  =  y(l)  =  =  3.2387945031585805075005735281910145 

Y(2)  =  y'(l)  =  Y(l)  cosh(l)  =  4.9977210779700644805224989588330930. 

Output  (2) :  With  a  =  0,  b  =  2  : 

Y(l)  =  y(2)  =  esinh(2)  =  37.5945994280492356058950044317377782 

Y(2)  =  y'(2)  =  Y(l)  cosh(2)  =  141.43823997622198995416647345901998. 

CASE2 :  ((y')2  +  l)y"'  -  3y'  (y")2  =  0;  y(l)  =  3,  f(l)  =  0,  y"(l)  =  -1/2; 

DE  =►  r(l)  =  0,  [4,  p.602, 7.11].  DY(3)  =  3Y(2)Y(3)2/[1  +  Y(2)2]. 

The  solution  is  given  by  y(t)  =  1  +  a/(3  —  t)  (1  +  t),  y'(t)  =  (1  —  t)/ (y  —  1), 

y"(t)  =  —  i  /  (y  —  l)  -  (y')V(y  -  !)•  yi  =  y2>  y 2  =  ys,  y3  =  f  =  3y2y3/(y2  + 1)- 

Output(l) :  With  a  =  1,  b  =  2  : 

Y(l)  =  yi(2)  =  y(2)  =  1  +  y/Z  =  2.7320508075688772935274463415058724 

Y(2)  =  yi(2)  =  f(2)  =  —1/^(3)  =  -.57735026918962576450914878050195746 

Y(3)  =  y"(2)  =  y"(2)  =  — 4/(3v/3)  =  -.76980035891950101934553170733594327. 

Output(2) :  With  a  =  1,  b  =  2.5  : 

Y(l)  =  yi(2.5)  =  y(2.5)  =  1  +  V7/2  =  2.3228756555322952952508078768196302 

Y(2)  =  yi(2.5)  =  f(2.5)  =  -3/V7  =  -1.1338934190276816816435496087025402 

Y(3)  =  y"  (2.5)  =  y"(2.5)  =  —32 /73^2  =  -1.7278375908993244672663613084991088. 

68.  (1264-1273) 

Call  DODE(F,  n,  Y,  T,  TOUT,  RERR,  AERR,  IND,  WK,  IWK) 

Double  precision  version  of  ODE  obtained  by  FCH.  Uses  same  arguments  except 
for  AERR=RERR=10 DPMPAR(l). 

69.  (1274— 1283)— Page 563 

Call  GERK(F,  n,  Y,  T,  TOUT,  RERR,  AERR,  IND,  GERROR,  WK,  IWK) 

ADAPTIVE  RKF  INITIAL  VALUE  DE  SOLVER. 

GERROR(i)  is  an  estimate  of  the  error  at  Y(i),  i  =  1, . . .  ,n.  The  remaining  variables 
play  the  same  role  as  in  ODE.  The  same  test  cases  used  in  ODE  are  used  with  GERK. 

70.  (1284-1293) 

Call  DGERK(F,  n,  Y,  T,  TOUT,  RERR,  AERR,  IND,  GERROR,  WK,  IWK) 

Double  precision  version  of  GERK  obtained  from  FCH.  Uses  same  arguments. 
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71.  ( 1 294— 1301) —Page  571 

Call  RK(n,  T,  h,  A,  F)— FOURTH  ORDER  RUNGE-KUTTA  . 

Solves  the  initial  value  problem  specified  in  ODE  by  using  the  Runge-Kutta  fourth  order 
procedure.  T  contains  the  initial  value  of  t  =  t0.  The  increment  by  which  t  is  increased 
at  each  call  of  RK  is  stored  in  h.  A(l), . . .  ,  A(n)  contain  the  values  yi(t0), . . .  ,yn(to)- 
A(n+  1),...  ,A(2n)  contain  the  values  y;  (t0),...  ,y£(to).  When  RK  is  called  the  values 
yi(t0  +  h), . . .  ,  yN(t0  +  h)  and  yi  (to  +  h), . . .  ,  y^(t0  +  h)  are  computed  and  stored  in 
A(l), . . .  ,  A(2n).  A(2n  +  1), . . .  ,  A(3n)  are  used  as  work  space.  F(t,  Z)  is  a  user  defined 
subroutine  which  computes  the  yj.  Two  CASES  are  carried  out. 

CASE1 :  9y2y"'  ~  45yy'y"  +  40  (f)3  =  0;  y(0)  =  1,  y*(0)  =  -3/2,  y"(0)  =  0; 

DE  =►  y"'(0)  =  15.  [4,  p.  602,  7.9]. 

Hence  for  yi(0)  =  A(l)  =  1,  y^O)  =  y2(0)  =  A(2)  =  —3/2, 

yf(0)  =  y'2(0)  =  A(3)  =  0,  Z(3)  =  [45Z(1)Z(2)Z(3)  -40Z(2)3]/[9Z(1)2]. 

The  solution  is:  y(t)  =  1/(1  + 1  +  5t2/4)3/2, 

?(t)  =  -(3/2)y5/3(l  +  St/2),  r(t)  =  5(y')2/(3y)  -  (15/4)ys/3- 
The  integration  is  carried  out  from  0  to  2  with  h=  1/32768,  where 
yi(2)  =  y(2)  =  (1/8 )3/2,  yi(2)  =  ?(2)  =  -9yi(2)/8,  yf(2)  =  y"(2)  =  105(l/8)3-5. 
Output :  A(l)  =  y(2)  =  .44194173824159220275052772631553065E(-1) 

A(2)  =  A  (4)  =  f( 2)  =  — ,49718445552179122809434369210497198E(— 1) 

A(3)  =  A(5)  =  y"(2)  =  .72506066430261220763758455098641747E(-1) 

CASE2 :  f  -  ty  =  0;  y(0)  =  1,  y'(O)  =  0,  DE  =»  y"(0)  =  0. 

The  solution,  derived  in  Appendix  B,  is  given  by 

7T  2 1  2 1^^ 

7(t) = wwm [V5 Ai(t) + Bi(t)1,  f (t) = s^ro/3)  7s  l2/3(_r : *•  (8) 

Ai  and  Bi  are  Airy  functions  [3,  p.446]  evaluated  from  Call  AI  and  Call  BI— Page 
(101-102).  I2/3 (t)  represents  a  Modified  Bessel  Function  which  can  be  evaluated  from 

Call  BESI— Page  97.  Z(2)=tZ(l). 

Output:  A(l)  =  y(2)  =  2.7308830178901459635915275691024881 
A(2)  =  A(3)  =  f( 2)  =  3.2595163616105247767976271320619734 

72.  (1302-1309) 

Call  DRK(n,T,h,  A,  F) -FOURTH  ORDER  RUNGE-KUTTA  . 

Double  precision  version  of  RK  obtained  from  FCH.  Uses  same  arguments. 

73.  (1310— 1317)— Page  573 
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Call  RK8(n,T,h,Y,DY,WK,F) -EIGHTH  ORDER  RUNGE-KUTTA  . 

Solves  the  initial  value  problems  specified  in  RK  by  using  the  eighth  order  Runge-Kutta 
procedure.  T  contains  the  initial  value  of  t  =  t0.  The  increment  by  which  t  is  increased 
at  each  call  of  RK8  is  denoted  by  h.  Y(l), . . .  ,  Y(n)  contain  the  values  yi(t0), .  • .  ,yn(to)- 
DY(1),...  ,  DY(n)  contain  the  values  y't  (t0), .. .  ,  y^(to)-  When  RK8  is  called  the  values 
yj(to  +  h), . . .  ,  yN(t0  +  h)  and  yi(t0  +  h), . . .  ,  y^(t0  +  h)  are  computed  and  stored  in 
Y(l), . . .  ,  Y(n),  DY(1), . . .  ,  DY(n).  WK  is  a  work  space,  and  F(t,  Z)  is  a  user  defined 
subroutine  which  computes  the  elements  of  DY  from  the  differential  equations.  The 
test  cases  are  the  same  as  those  used  in  RK  above. 

74.  (1318-1325) 

Call  DRK8(n,  T,  h,  Y,  DY,  WK,  F)-EIGHTH  ORDER  RUNGE-KUTTA  . 

Double  precision  version  of  RK8  obtained  from  FCH.  Uses  same  arguments. 
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III.  Testing  Of  Bit  Operations 

The  following  intrinsic  bit  functions  were  tested : 

1.  BTEST  2.  IBITS  3.  IBCLR  4.  IBSET  5.  ISHFTC 

6.  ISHFT  7.  IAND  8.  MVBITS  9.  IOR  10.  IEOR 

In  case  an  error  occurs,  the  error  message  will  be  listed  with  an  identifying  symbol 
<  KX  >,  where  K  denotes  an  integer  and  X  a  letter.  The  symbol  will  also  be  included  in 
a  comment  statement  of  the  code  to  help  locate  the  section  of  the  code  that  was  operating 
when  the  error  occured. 

Let  MSIZE  denote  the  number  of  bits  in  an  integer  word.  Let  ISIZE  =  min(MSIZE,  128). 
If  word  size  exceeds  128  bits  only  the  first  128  bits  are  tested.  Let  I  be  an  integer  with  ISIZE 
alternating  zero  and  one  bits  with  bit  zero  set  to  one,  and  let  J  denote  its  complement. 
Tests  2,  3,  and  4  are  only  done  if  ISIZE  is  even. 

1. The  following  test  is  done  twice— once  with  I  and  once  with  J. 

For  each  bit  position  K  (K=  0, . . .  ,  ISIZE), 

let  N  be  the  value  of  the  bit  at  K,  and  let  M  be  the  value  of  the  bit  at  K  detected  by 
BTEST. 

Verify  that  M  =  N  and  that  I  is  unchanged.  <  1A  > 

Use  IBCLR  to  dear  the  bit  at  K  in  I.  Use  IBITS  to  check  that  it  is  cleared  and  I 
is  unchanged. <  IB  > 

Also  use  IBITS  to  ascertain  that  the  result  is  unchanged  except  at  K.<  IB  > 

Test  IBSET  as  was  done  with  IBCLR.  <  1C  > 

2.  For  each  L(L  =  0, 1, 2, . . .  ,  ISIZE) 

If  L  is  odd  let  M  =  J 
if  L  is  even  let  M  =  I 

for  each  II  (II  =  L,  L  +  1, . . .  ,  ISIZE)  and  II  ^  0 
let  K  =  ISHFTC(I,  L,  II)  (left  circular  shift) 
compute  K  using  BTEST,  IBCLR,  IBSET 
if  the  two  results  differ  print  an  error  message<  2A  > 

let  K  =  ISHFTC(I,  — L,  II)  (right  circular  shift)  and  repeat  above  tests<  2B  > 

let  K  =  ISHFTC  (M,  1,  ISIZE)  and  let  MM  =  complement  of  M 
K  ^  MM  print  error  message  <  2C  > 

let  K  =ISHFTC(M,-1, ISIZE) 
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if  K  7^  MM  print  error  message  <  2D  > 

3.  Let  (McJ)  denote  the  M  corresponding  bits  of  J 
let  KI  =  I  and  let  KJ  =  J 

for  each  L  (L  =  0, 1, . . .  ,  ISIZE  —  1) 
if  L  ^  0  then  use  IBCLR  to: 

clear  the  right  most  L  bits  of  KI 
clear  the  left  most  L  bits  of  KJ 

ISIZE— L  L 

K  =  ISHFT(I1,  L)  =  IT •  ^ ■  lll'bO  •  -  00(J  II  contains  all  one  bits 

((ISIZE— L)cl>  l 

K  =  IAND(K,  I)  =  xxx  -x^OO^OOO  =  KI 
if  K  ^  KI  print  error  message<  3A  > 

L  ISIZE— L 

K  =  ISHFT(I1,  -L)  =  ^O^OO^TT^TTT 

L  ((ISIZE— L)cJ) 

K  =  IAND(K,  J)  =  b0- -000' xxx  — xxx  =  KJ 
if  K  7^  KJ  print  error  message<  3B  > 

4.  For  each  L  (L  =  0, 1, . . .  ,  ISIZE  —  1 ) 

For  each  N  (N  =  0, . . .  ,  ISIZE  -  L) 

K  =  J 

N  bits  are  moved  from  I  to  K  beginning  with  bit  position  L  in  both  I  and  K 
call  MVBITS(I,L,N,K,L) 

(McJ)  (Ncl)  (LcJ) 

K  =  xx  •  ••  xxx' xx  •  —  xxx  xx  •  —  xxx,  M  +  N  +  L  =  ISIZE 


KI  =  IEOR(K,  I),  KJ  =  IEOR(K,  J) 

M  N  L 

KI  =  ll  — llfOO  — 000'll  —  lll' 

M  N  L 

KJ  =  {xT^OOCTl  1  —  1 1  foo  •  -  000' 

if  IEOR(KI,  KJ)  7^  II  (all  ones)  print  an  error  message<  4A  > 
if  IAND(KI,  KJ)  7^  0  print  an  error  message<  4B  > 


KI  =  IOR(K,  I),  KJ  =  IOR(K,  J) 

(Ncl)  L 


KI  =  11  *  -  *  111: 


■xxxll  •  •  •  111 
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(McJ)  N  (LcJ) 

kj  =  sor^^lT^nTsar^ioS 

if  IAND(KI,  KJ)  ^  K  print  an  error  message<  4C  > 


IV.  Testing  Of  String  Operations 

The  following  intrinsic  string  functions  were  tested : 

1.  ACHAR  2.  CHAR  3.  IACHAR  4.  ICHAR  5.  ADJUSTL 

6.  ADJUSTR  7.  LEN  8.  LEN_TRIM  9.  INDEX  10.  LLT 

11.  LLE  12.  LGE  13.  LGT  14.  REPEAT  15.  TRIM 

16.  SCAN  17.  VERIFY 

In  case  an  error  occurs,  the  error  message  will  be  listed  with  an  identifying  symbol 
<  KX  >,  where  K  denotes  an  integer  and  X  a  letter.  The  symbol  will  also  be  included  in 
a  comment  statement  of  the  code  to  help  locate  the  section  of  the  code  that  was  operating 
when  the  error  occurred. 

Let  AA  = ’abc- •  •  xyz’  BB  =  ’ABC---  XYZ’ 

CC  =  '1234567890'—  =  \~ !  @  #  $  %  *  &  *()  - 
ALL  (1+ 1:1+1  )=CHAR(I)  I  =  0, 1, . . .  ,  127 

1.  For  1  =  0,1,...  ,127 

C  =  ACHAR(I)  D  =  CHAR(I) 
if  C  7^  D  print  error  message  <  1A  > 

J  =  IACHAR(D)  K  =  ICHAR(D) 
if  J  ^  I  print  error  message<  IB  > 
if  K  7^  I  print  error  message  <  1C  > 

2.  ABC  =  AA//BB//CC 

for  I  =  0, 1, . . .  ,72 
if  I  <  72  then 
ABCX  =  ABCY  =  ABC 
if  I  >  1  ABCX(1  : 1)  ='  b'  b  s  a  blank 
ABCX  =  ADJUSTL(ABCX) 

ABCX  has  I  blanks  on  the  right;  its  length  (minus  trailing  blanks)  is  72  - 1. 
J=LEN_TRIM(ABCX)  K=LEN(ABCX(1:J)) 
if  J  ^  72  —  I  print  error  message<  2A  > 
if  J  7^  K  print  error  message<  2B  > 
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in  ABCY  move  character  positions  1+1,. . .  ,72  into  character 
positions  1,2,...  ,72  —  1,  then  blank  out  positions  72  —  I  +  1, . . .  ,72 
if  ABCX  /  ABCY  print  error  message<  2C  > 

if  I  >  1 ,  then 
ABCX  =  ABCY  =  ABC 
ABCX(1 :  72)  ='  b' 

ABCX  =  ADJUSTR(ABCX) 

in  ABCY  move  character  positions  1,2,...  ,1  —  1  to  character 
positions  74  —  I, . . .  ,72  and  blank  out  positions  1, 2, . . . ,  73  —  I 
if  (ABCX  7^  ABCY)  print  error  message<  2D  > 

3.  DD  =  CC  DD(1:1)  =  V  ABC  =  DD//DD//DD  DD(2:2)  =  V 

the  arguments  of  INDEX  are  INDEX(string,  substring,  back).  If  back  is  missing  it 

is  equivalent  to  back  =  .FALSE. 

when  string  is  shorter  than  substring  INDEX  =  0 

J  =  INDEX(DD(2  :  72),  ABC) 

K  =  INDEX(DD(2  :  72),  ABC,  .FALSE.)  N  =  INDEX(DD(2  :  72),  ABC,  .TRUE.) 
if  J  ^  K  print  error  message  <  3 A  > 
if  J  0  print  error  message<  3B  > 
if  N  7^  0  print  error  message<  3C  > 

For  each  I  (I  =  0, 1 , . . . ,  71 ) 

For  each  II  (II  =  1, 2, . . . ,  I  +  1) 

J  =  INDEX(ABC,  DD(II:I)) 

K  =  INDEX  (ABC,  DD(II:I),  .FALSE.) 

N  =  INDEX(ABC,  DD(II:I)),  .TRUE.) 
if  J  K  print  error  message  <  3D  > 

if  II  >  I  then 

DD(II:I)  is  empty,  hence 

if  back  =  .FALSE,  then  INDEX=  1 

if  back  =  .TRUE.  INDEX  =  73. 

therefore  if  J  ^  1  print  error  message<  3E  > 

if  N  7^  73  print  error  message  <  3F  > 

if  II  =  1  and  I  >  1  then 
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DD(II:I)  begins  with  2  blanks,  hence  DD(II:I)  is  not  in  ABC,  INDEX  =  0 
if  J  ^  0  print  error  message<  3G  > 
if  N  ^  0;  print  error  message<  3H  > 

if  II  =  2  and  1^2  then 

DD(II:I)  begins  with  'b3',  this  substring  doesn’t  occur  in  ABC,  since  in 
ABC  '3'  is  always  preceded  by  '2' 
if  J  ^  0  print  error  message<  31  > 
if  N  ^  0  print  error  message<  3J  > 

else 

II  <  I  and  (II  >  1  or  I  <  1)  and  (11^2  or  1=2) 
therefore  II  =  I  =  1  or  II  =  I  =  2  or  2  <  II  <  I 

if  II  =  2  then  DD(II :  I)  ='  b'  occurs  at  ABC(1:1),  ABC(27:27),  ABC(53:53) 
print  error  message  if  J  ^  1  <  3K  > 
print  error  message  if  N  ^  53  <  3L  > 

if  II  ^  2  then  DD(II:I)  occurs  at  ABC(M:M),  where  M  =  II,  II  +  26,  II  +  52 
print  error  message  if  J  ^  II  <  3K  > 
print  error  message  if  N  ^  II  +  52.  <  3L  > 

4.  For  each  J  (J  =  1, 2, . . .  ,  128) 

For  each  I  (1=1,2,...,  128) 

let  II  =  ALL(I :  I)  let  JJ  =  ALL(J  :  J) 
if  LLT(II,  JJ)  is  TRUE  and  I  >  J  print  error  message<  4A  > 
if  LLE(II,  JJ)  is  TRUE  and  I  >  J  print  error  message<  4B  > 
if  LGE(II,  JJ)  is  TRUE  and  I  <  J  print  error  message<  4C  > 
if  LGT(II,  JJ)  is  TRUE  and  I  <  J  print  error  message<  4D  > 

5.  For  each  I  (I  =  1, 2, . . .  ,  8) 

For  each  K  (K  =  0, 1, . . . ,  8) 

DD  =  AA  set  last  K  characters  of  DD  to  blank  M  =  26  —  K 
if  TRIM(DD)  ^  DD(1  :  M)  print  error  message<  5 A  > 
if  LEN(TRIM(DD))  ^  M)  print  error  message<  5A  > 

if  K  >  0  then 

let  ABCX  =  AA(1  : 1)//AA(1  : 1)  •  •  •  K  times 
let  ABCY  =  REPEAT(AA(1  : 1),  K) 
if  ABCX  ^  ABCY  print  error  message<  5B  > 
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6.  ABC  =  AA//BB 

For  each  I  (I  =  1,2,...  ,52) 

K  =  SCAN(ABC(1: 52),  AA) 

M  =  SCAN (ABC (1 : 52),  AA,  .FALSE.) 

N  =  SCAN (ABC (1 : 52),  AA,  .TRUE.) 

M(N)  =  the  position  number  of  the  left  most  (right  most)  character 
in  ABC(1 : 52)  that  is  in  AA 
if  K  ^  M  print  error  message<  6A  > 

if  I  <  26  then 

if  M  7^  1  print  error  message  <  6B  > 

if  N  7^  26  —  (I  —  1)  print  error  message<  6C  > 

if  I  >  26  then 

if  M  7^  0  print  error  message<  6D  > 
if  N  ^  0  print  error  message<  6E  > 

K  =  VERIFY (ABC(1 : 52),  AA) 

M  =  VERIFY  (ABC  (1 : 52),  AA,  .FALSE.) 

N  =  VERIFY (ABC(1 : 52),  AA,  .TRUE.) 

M(N)  =  the  position  number  of  the  left  most  (right  most)  character 
in  ABC(1 : 52)  that  is  not  in  AA 
if  K  ^  M  print  error  message<  6F  > 

if  N  ^  52  —  (I  —  1)  print  error  message  <  6G  > 
if  I  <  26  then 

if  M  7^  27  —  (I  —  1)  print  error  message<  6H  > 
if  I  >  26  then 

if  M  7^  1  print  error  message<  61  > 
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V.  How  To  Use  TEST90 

TEST90  consists  of  four  Fortran  90  source  files— TEST90.F,  LIB.F,  LIBS.F,  and  LIBD.F. 
TEST90.F  contains  the  driver  program,  the  bit  and  character  test  routines,  the  block  data 
routines  containing  the  safety  factors  (SF)  and  stored  epsilons  (SE)4,  and  some  miscellaneous 
routines.  It  was  successfully  tested  on  Fortran  90  compilers  of  the  CRAY  EL98  and  the  IBM 
PC. 

LIB.F  contains  74  basic  testing  routines  and  some  supporting  software.  Each  is  used  to 
test  one  of  the  74  subprograms  in  Section  II.  The  supporting  routines  to  LIB.F  are  in  LIBS.F 
and  LIBD.F.  The  routines  tested  in  LIBS.F  are  all  MLIB  routines.  The  routines  tested  in 
LIBD.F  are  all  double  precision  routines  created  by  FCH. 

The  first  routine  in  LIBS.F  is  IPMPAR.  The  machine  constants  in  IPMPAR  must  be 
set  for  the  particular  compiler  on  the  particular  machine  being  used.  Instructions  are  given 
in  the  in-line  documentation  of  IPMPAR  for  defining  the  constants  that  are  needed.  If 
constants  are  not  provided  for  the  machine  being  used  then  consult  the  Fortran  compiler 
manual  for  that  machine  or  the  appendix  of  [1].  These  constants  in  LIBS.F  were  originally 
set  for  the  64  bit  CRAY  EL98.5 

To  use  TEST90  one  needs  to  compile  and  link  the  four  Fortran  90  source  files-TEST90.F, 
LIB.F,  LIBS.F,  and  LIBD.F  and  run  the  resulting  executable  code. 

If  a  routine  of  LIB.F  reveals  an  error,  a  print  out  of  the  case  involved  gives  the  true  and 
calculated  values,  the  relative  error  EPS  or  EPD,  and  the  case  number  N(1  <  N  <  1325).  If 
a  routine  in  LIB.F  encounters  no  errors,  then  ’x  ok’  is  printed  where  x  is  the  routine  name. 

If  a  bit  or  sting  error  occurs,  a  print  out  gives  the  intrinsic  function  with  the  input,  the 
true  and  calculated  values,  and  the  error  message  symbol  <  KX  >  6.  If  no  bit  errors  occur 
then  ’bit  functions  ok’  is  printed.  If  no  string  errors  occur  then  ’string  functions  ok’  is 
printed.  If  no  errors  are  discovered  using  LIB.F  and  all  the  bit  function  and  string  function 
tests,  then  ’everything  ok’  is  printed. 


4  Also  listed  in  Appendix  C 

5 Since  none  of  the  routines  in  MLIB  use  integers  that  exceed  46  bits  either  compiler  option  (46  or  64  bits)  may  be  specified 

when  using  the  CRAY  EL98. 

6  See  Sections  III  and  IV. 
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VI.  Index 

An  integer  without  parentheses  in  this  index  gives  the  number  of  the  page  where  a  callable 
routine  is  discussed.  An  integer  in  parentheses  refers  to  the  corresponding  page  number  in 
[1].  The  seventeen  double  precision  directly  callable  routines  obtained  by  using  FCH  are 
identified  by  a  superscript  *.  The  eleven  callable  routines  used  as  supporting  routines  are 
marked  by  the  superscript  **. 


ARTNQ  (11),  5 
BESI**  (97),  25 
BRATIO  (88),  8 
BSSLJ  (92),  10 
BSSLK  (99),  11 
CBESK  (99),  12 
CBRT  (11),  5 
CBSPL  (455),  22 
CBSSLI**  (96),  12 
CBSSLJ  (92),  11 
CBSSLK**  (98),  12 
CBSSLY**  (95),  26 
CERF  (45),  7 
CERFC  (46),  7 
CHEBY  (433),  20 
CIRCV**  (125),  24,  25 
CSINT  (465),  24 
DARTNQ  (11),  5 
DAXPY  (175),  15 
DBESK  (98),  12 
DBRAT  (89),  10 
DBSLJ*  (92),  11 
DBSSLI**  (96),  13 
DBSSLJ  (91),  11 
DBSSLK*  (99),  12,  13 
DBSSLY**  (94),  26 


DCBRT  (11),  5 
DCBSPL*  (455),  24 
DCERF  (47),  7 
DCERFC  (48),  8 
DCPOLY  (157),  14 
DCSINT*  (465),  24 
DCSQRT  (13),  5 
DEIGV  (357),  17 
DELLPI**  (109),  20 
DERF  (48),  7 
DERFC1  (49),  7 
DGERK*  (563),  27 
DHTRP*  (447),  21 
DKROUT*  (216),  15 
DLLSQ  (385),  18 
DMSLV  (217),  15 
DODE*  (551),  27 
DOPTF*  (409),  18 
DPFIT*  (451),  22 
DQAGI  (540),  26 
DQAGS  (534),  25 
DQSUBA*  (533),  26 
DREXP  (29),  6 
DRK*  (571),  28 
DRK8*  (573),  29 
DRLOG  (32),  6 


DRPOLY  (157),  13 
DSFVAL*  (525),  25 
DSLVMP*  (221),  16 
DSMPLX*  (413),  19 
DSSPLX*  (413),  20 
DZERO  (151),  13 
EIGV  (353),  16 
ERF  (46),  6 
ERFC1  (46),  7 
GERK  (563),  27 
GRA.DL**  (524),  24 
HTRP  (447),  20 
KROUT  (216),  15 
LLSQ  (377),  17 
MSLV  (217),  15 
ODE  (551),  26 
OPTF  (409),  18 
PCOEFF**  (449),  21 
PFIT  (451),  21 
QAGI  (539),  26 
QAGS  (531),  25 
QSUBA  (533),  26 
REXP  (29),  6 
RK  (571),  28 
RK8  (573),  29 
RLOG  (31),  6 
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SAXPY  (175),  14 
SFVAL  (525),  24 
SLVMP  (221),  15 


SMPLX  (413),  18  ZEROIN  (151),  12 

SSPLX  (413),  20 
TRMESH**  (523),  24 
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VII.  Summary 

The  objective  of  this  work  was  to  use  the  NSWC  Mathematical  Library  (MLIB)[1]  as  an 
aid  in  testing  a  new  Fortran  90  compiler  used  with  the  CRAY  EL98  mainframe  computer. 
One  can  also  view  the  testing  as  an  aid  in  establishing  whether  MLIB  can  operate  in  a 
Fortran  90  environment.  The  final  result  is  a  transportable  driver  program,  TEST90,  which 
is  not  restricted  to  testing  only  the  new  90-compiler.  Any  post  Fortran  77  compiler  can 
be  tested  provided  it  is  limited  to  a  single  precision  word  length  of  64  bits  or  less.  This 
limitation  is  imposed  because  the  routines  in  MLIB  are  based  on  a  64  bit  word  length  with 
a  48  bit  floating  point  mantissa. 

TEST90  has  been  used  to  test  the  latest  Fortran  90  compilers  of  the  CRAY  EL98  and 
IBM  PC  machines;  no  errors  were  found,  however  TEST90  did  reveal  a  complex  arthmetic 
error  in  an  earlier  version  of  the  Cray  EL98  compiler.  MLIB  routines  ran  under  TEST90 
without  any  problems  on  both  machines. 
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APPENDIX  A 

INVERSE  OF  THE  HILBERT  MATRIX  OF  ORDER  10 


This  appendix  contains  the  elements  of  the  inverse  of  the  tenth  order  Hilbert  matrix.  The 
matrix  is  given  by  rows.  The  jth  row  is  ended  with  ■  (j). 


100,  -4950,  79200,  -600600,  2522520,  -6306300,  9609600,  -8751600,  4375800,  -9237801  (1) 

-4950,  326700,  -5880600,  47567520,  -208107900, 

535134600,  -832431600,  770140800,  -389883780,  831402001  (2) 

79200,  -5880600,  112907520,  -951350400,  4281076800, 

-11237826600,  17758540800,  -16635041280,  8506555200,  -18290844001(3) 

-600600,  47567520,  -951350400,  8245036800,  -37875637800, 

101001700800,  -161602721280,  152907955200,  -78843164400  ,170714544001  (4) 

2522520,  -208107900,  4281076800,  -37875637800,  176752976400, 

-477233036280,  771285715200,  -735869534400,  382086104400,  -832233402001(5) 

-6306300,  535134600,  -11237826600,  101001700800,  -477233036280, 

1301544644400,  -2121035716800,  2037792556800,  -1064382719400,  2330253525601(6) 

9609600,  -832431600,  17758540800,  -161602721280,  771285715200, 

-2121035716800,  3480673996800,  -3363975014400,  1766086882560,  -3883755876001  (7) 

-8751600,  770140800,  -16635041280,  152907955200,  -735869534400, 

2037792556800,  -3363975014400,  3267861442560,  -1723286307600,  3804495552001  (8) 

4375800,  -389883780,  8506555200,  -78843164400,  382086104400, 

-1064382719400,  1766086882560,  -1723286307600,  912328045200,  -2021138262001  (9) 
-923780,  83140200,  -1829084400,  17071454400,  -83223340200, 

233025352560,  -388375587600,  380449555200,  -202113826200,  449141836001  (10) 
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APPENDIX  B 

DERIVATION  OF  EQUATION  (8) 


The  objective  here  is  to  derive  the  solution  of  the  Airy  differential  equation  (DE)  with 
specified  initial  values,  where 

At)  - 1  y(t)  =  0;  y(0)  =  1 y'(0)  =  0.  (B-l) 

Linearly  independent  solutions  of  the  DE  in  (B-l)  or  (8)  are  the  Airy  functions  [Bl,  p.446] 
denoted  by  Ai(t)  and  Bi(t).  Thus 

y  =  cj  Ai(t)  +  c2  Bi(t). 


The  initial  value  conditions  yield 

ci  Ai(0)  +  c2  Bi(0)  =  1,  ci  Ai'(0)  +  c2  Bi'(0)  =  0.  (B-2) 

The  solution  of  (B-2)  is  given  by 

d  =  Bi'(0)/W,  c2  =  —  Ai'(0)/W, 
where  W  denotes  the  Wronskian  of  (B-l), 

W  =  Ai(0)  Bi'(0)  —  Ai'(0)  Bi(0). 

From  [Bl,  p.446],  we  have 


Ai(0)  =  Bi(0)/V3  =  3-2/3/r(2/3),  Ai'(0)  =  -Bi'(0)/V3  =  -3-1/3/r(l/3). 

Therefore  using  6.1.17  of  [Bl,  p.256]  W  =  1/tt;  y(t),  y'(t),  after  some  algebra,  are  given  by 

y(t)  =  3V3r(i/3)[VSAi(t)  +  Bi(t)1’  m  =  3V3r(i/3)lv^Ai'(t)  +  Bi'(t)]-  (B_3) 

The  quantities  Ai'(t)  and  Bi'(t)  can  be  expressed  in  terms  of  Modified  Bessel  functions 
[Bl,  p.447]  as 


Ai(t)  =  3 


,2t3/2 


2 13/2 , 


I2/3  (  o  )  ~  I — 2/3  (  o  ) 


•  Bi'(t)  =  7!  L 


,2t3/2 


2t3/2, 


I2/3  (  q  )  +  1-2/3  (  o  ) 


which  gives,  using  (B-3), 

,  7T  2t  _  ,2t3/2 

y^)- 31/3  r(i/3)V3 12/3  3 

Bl-Abramowitz,  M.,  Stegun,  I.,  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS,  Dover  Publications, N.Y.,  1965. 
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APPENDIX  C 

RELATIVE  ERROR  TYPE (7^.  or  T ),  SE  AND  SF  FACTORS 

This  appendix  contains  for  each  case  of  Section  II :  its  number,  the  relative  error  type 
(either  71  or  T ),  the  SE  factor,  and  the  SF  factor  (see  (3)  and  (5)  of  Section  I). 


1 

n 

4.75E01 

5.00E00 

2 

n 

1.22E01 

3.33E00 

3 

71 

1.25E00 

5.00E00 

4 

n 

1.00E00 

5.00E00 

5 

n 

1.25E00 

5.00E00 

6 

71 

1.16E01 

3.29E00 

7 

71 

1.25E00 

3.63E01 

8 

n 

8.61  E01 

3.00E00 

9 

71 

1.03E01 

4.07E00 

10 

n 

1.25E00 

5.00E00 

11 

n 

1.00E00 

5.00E00 

12 

71 

1.25E00 

5.00E00 

13 

n 

1.09E01 

3.56E00 

14 

71 

2.13E01 

3.00E00 

15 

71 

1.25E00 

5.00E00 

16 

n 

1.25E00 

5.00E00 

17 

71 

1.25E00 

5.00E00 

18 

71 

1.25E00 

5.00E00 

19 

n 

1.25E00 

5.00E00 

20 

71 

1.25E00 

5.00E00 

21 

71 

1.25E00 

5.00E00 

22 

n 

1.25E00 

5.00E00 

23 

71 

1.25E00 

5.00E00 

24 

71 

1.25E00 

5.00E00 

25 

n 

1.25E00 

5.00E00 

26 

71 

1.25E00 

5.00E00 

27 

71 

1.25E00 

5.00E00 

28 

n 

1.25E00 

5.00E00 

29 

71 

1.25E00 

5.00E00 

30 

71 

1.25E00 

5.00E00 

31 

n 

1.25E00 

5.00E00 

32 

71 

1.25E00 

5.00E00 

33 

71 

1.25E00 

5.00E00 

34 

n 

1.25E00 

5.00E00 

35 

71 

1.25E00 

5.00E00 

36 

71 

1.25E00 

5.00E00 

37 

n 

1.25E00 

5.00E00 

38 

71 

1.25E00 

5.00E00 

39 

71 

1.25E00 

5.00E00 

40 

n 

1.25E00 

5.00E00 

41 

71 

1.25E00 

5.00E00 

42 

71 

1.25E00 

5.00E00 

43 

n 

1.00E00 

5.00E00 

44 

71 

1.46E02 

3.00E00 

45 

71 

1.25E00 

5.00E00 

46 

n 

1.25E00 

5.00E00 

47 

71 

1.25E00 

5.00E00 

48 

71 

1.25E00 

5.00E00 

49 

n 

1.00E00 

5.00E00 

50 

71 

6.43  E01 

1.49E01 

51 

71 

1.25E00 

5.00E00 

52 

n 

1.25E00 

5.00E00 

53 

71 

1.25E00 

5.00E00 

54 

71 

1.25E00 

5.00E00 

55 

n 

1.25E00 

5.00E00 

56 

71 

1.25E00 

5.00E00 

57 

71 

1.25E00 

5.00E00 

58 

n 

1.25E00 

5.00E00 

59 

71 

1.25E00 

5.00E00 

60 

71 

1.25E00 

5.00E00 

61 

n 

1.25E00 

5.00E00 

62 

71 

1.25E00 

5.00E00 

63 

71 

1.25E00 

5.00E00 

64 

n 

1.25E00 

5.00E00 

65 

71 

1.25E00 

5.00E00 

66 

71 

1.25E00 

5.00E00 

67 

n 

1.25E00 

5.00E00 

68 

71 

1.25E00 

5.00E00 

69 

71 

1.25E00 

5.00E00 

70 

n 

1.25E00 

5.00E00 

71 

71 

1.25E00 

5.00E00 

72 

71 

1.25E00 

5.00E00 

73 

n 

1.25E00 

5.00E00 

74 

71 

1.25E00 

5.00E00 

75 

71 

1.25E00 

5.00E00 

76 

n 

1.25E00 

5.00E00 

77 

71 

1.25E00 

5.00E00 

78 

71 

1.25E00 

5.00E00 

79 

n 

1.25E00 

5.00E00 

80 

71 

1.25E00 

5.00E00 

81 

71 

1.25E00 

5.00E00 

C-3 


NSWCDD/TR-98/75 


82  72.  1.25E00  5.00E00 
85  n  1.25E00  5.00E00 
88  72. 1.11E01  5.00E00 
91  72.  1.25E00  5.00E00 
94  72.  2.60E01  3.33E00 
97  72  1.25E00  5.00E00 
100  ft  1.25E00  5.00E00 
103  72  1.25E00  5.00E00 
106  72  1.25E00  5.00E00 
109  72  1.25E00  5.00E00 
112  72. 4.99E01  5.00E00 
115  72  1.25E00  5.00E00 
118  72  1.25E00  5.00E00 
121  72  1.25E00  5.00E00 
124  72  1.75E02  3.00E00 
127  72  1.25E00  4.44E01 
130  72  1.25E00  5.00E00 
133  72  1.25E00  5.00E00 
136  72  3.97E01  3.00E00 
139  ft  1.25E00  5.00E00 
142  ft  5.63E00  5.00E00 
145  ft  1.25E00  5.00E00 
148  ft  1.02E01  5.00E00 
151  ft  1.25E00  5.00E00 
154  ft  1.25E00  5.00E00 
157  ft  1.25E00  5.00E00 
160  ft  4.96E02  5.00E00 
163  ft  1.25E00  5.00E00 
166  ft  1.25E00  5.00E00 
169  ft  1.46E01  4.80E00 
172  ft  1.78E02  3.00E00 
175  ft  1.06E01  3.18E00 
178  ft  1.25E00  5.00E00 
181  ft  1.25E00  5.00E00 


83  ft  1.25E00  5.00E00 
86  ft  3.24E01  3.00E00 
89  ft  1.25E00  5.00E00 
92  ft  1.25E00  4.35E01 
95  ft  1.25E00  5.O0E00 
98  ft  1.25E00  5.00E00 
101  ft  1.25E00  5.00E00 
104  ft  1.25E00  5.00E00 
107  ft  1.25E00  5.00E00 
110  ft  1.25E00  5.00E00 
113  72. 1.25E00  5.OOEO0 
116  72. 5.12E00  5.00E00 
119  72-  1.25E00  5.00E00 
122  ft  9.33E00  5.00E00 
125  72. 1.25E00  5.00E00 
128  72. 4.75E02  5.00E00 
131  ft  9.60E00  5.00E00 
134  ft  1.25E00  5.00E00 
137  ft  1.25E00  5.00E00 
140  ft  1.61E02  3.00E00 
143  ft  1.25E00  5.00E00 
146  ft  7.00E00  5.00E00 
149  ft  7.98E00  5.00E00 
152  ft  3.90E01  3.00E00 
155  ft  1.25E00  5.00E00 
158  ft  1.80E01  3.00E00 
161  ft  1.25E00  5.00E00 
164  ft  1.29E01  3.00E00 
167  ft  1.25EOO  5.00E00 
170  ft  1.25E00  5.00E00 
173  72  9.17E00  8.41E00 
176  72. 4.91E02  5.00E00 
179  ft  1.25E00  8.03E01 
182  ft  1.25E00  5.00E00 


84  ft  5.94E00  5.00E00 
87  ft  1.25E00  5.00E00 
90  ft  1.25E00  5.00E00 
93  ft  1.25E00  5.00E00 
96  ft  1.90E01  4.04E00 
99  ft  1.25E00  5.00E00 
102  ft  1.25E00  5.00E00 
105  ft  1.25E00  5.00E00 
108  ft  1.25E00  5.00E00 
111  72. 1.25E00  5.00E00 
114  ft  1.25E00  5.00E00 
117  72. 7.10E00  5.00E00 
120  ft  3.92E01  3.00E00 
123  ft  7.83E00  5.00E00 
126  ft  1.25E00  2.17E01 
129  ft  1.25E00  5.00E00 
132  ft  5.87E00  5.00E00 
135  ft  1.25E00  5.00E00 
138  ft  1.25E00  5.00E00 
141  ft  1.25E00  5.00E00 
144  ft  5.00E02  5.00E00 
147  ft  1.25E00  5.00E00 
150  ft  1.25E00  5.00E00 
153  ft  1.25E00  5.00E00 
156  ft  1.71E02  3.00E00 
159  ft  1.62E01  3.00E00 
162  ft  6.36E00  5.00E00 
165  ft  1.25E00  2.50E01 
168  ft  3.94E01  3.00E00 
171  ft  7.38E00  5.00E00 
174  ft  1.43E01  3.41E00 
177  ft  1.25E00  6.32E01 
180  ft  7.45 E01  3.00E00 
183  ft  6.55E00  5.00E00 
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184  TZ  1.25E00  5.00E00 
187  TZ  1.25E00  5.00E00 
190  TZ  1.11E01  5.00E00 
193  72  1.25E00  5.00E00 
196  72  2.74E02  3.00E00 
199  72  2.51E01  3.00E00 
202  72  6.32E00  5.00E00 
205  72  1.25E00  5.00E00 
208  72  7.19E01  5.00E00 
211  TZ  1.25E00  3.16E01 
214  72  9.01EOO  3.00E00 
217  72  9.52E00  5.00E00 
220  72  6.79E00  5.00E00 
223  72  1.25E00  5.00E00 
226  TZ  1.25E00  5.00E00 
229  72  1.25E00  5.00E00 
232  72  8.70E02  3.00E00 
235  72  1.39E01  5.00E00 
238  TZ  1.25E00  5.00E00 
241  72  1.46E01  5.00E00 
244  72  2.48E02  3.00E00 
247  72  1.25E00  1.25E02 
250  72  6.82E00  5.62E00 
253  72  8.21  EOO  5.00E00 
256  72  1.25 EOO  5.00E00 
259  72  1.00 EOO  5.00E00 
262  72  l.OOEOO  5.00E00 
265  72  7.84E00  5.00E00 
268  72  1.25E00  5.00E00 
271  72  1.25 EOO  5.00E00 
274  72  1.25E00  5.00E00 
277  TZ  l.OOEOO  5.00 EOO 
280  72  l.OOEOO  5.00E00 
283  7 2 1.01E01  5.00E00 


185  72  1.25 EOO  5.00E00 
188  72  7.32E00  5.00E00 
191  72  1.25 EOO  5.00E0O 
194  TZ  1.36E01  5.00E0O 
197  72  6.13E00  5.00E00 
200  72  8.55E02  3.00E00 
203  72  1.01E01  5.00E00 
206  72  1.25E00  5.00E00 
209  72  1.25E00  5.00E00 
212  72  2.65E02  3.00EOO 
215  U  3.10E01  5.00 EOO 
218  n  1.87E01  5.00E00 
221  n  1.29E01  5.00E00 
224  n  7.41  E01  3.00E00 
227  Tl  1.25E00  5.00E00 
230  U  1.25E00  2.40E01 
233  n  1.25 EOO  2.76E01 
236  TZ  2.07E01  5.00E00 
239  U  1.25E00  5.00E00 
242  Tl  1.25E00  5.00E00 
245  1Z  1.58E01  5.00E00 
248  TZ  8.08E02  3.00E00 
251  Tl  1.02E01  3.09E00 
254  TZ  1.25 EOO  5.00E00 
257  TZ  1.25 EOO  5.00E00 
260  TZ  1.25 EOO  5.00E00 
263  TZ  1.25 EOO  5.00E00 
266  TZ  1.25 EOO  5.00E00 
269  TZ  1.25E00  5.00E00 
272  TZ  1.25E00  2.39E01 
275  TZ  1.25E00  5.00E0O 
278  TZ  1.25 EOO  5.00E00 
281  TZ  1.25 EOO  5.00E00 
284  TZ  5.08E00  5.00E00 


186  TZ  1.25E00  5.00E00 
189  72. 9.76E00  5.00E00 
192  72. 7.09E01  3.00E00 
195  TZ  5.22E00  5.00E00 
198  TZ  1.77E01  3.00E00 
201  TZ  1.25E00  5.00E00 
204  72  7.33E00  5.00E00 
207  72  1.25E00  5.00E00 
210  72  5.07E00  5.00E00 
213  72  1.25E00  2.16E01 
216  72  8.12E02  3.00E00 
219  72  6.18E00  5.00E00 
222  72  1.25E00  5.00E00 
225  72  1.25E00  5.00E00 
228  72  2.76E02  5.00E00 
231  72  1.25E00  1.26E02 
234  72  7.94E00  5.00E00 
237  72  1.49E01  4.47E00 
240  72  6.72E01  5.00E00 
243  72  1.25E00  5.00E00 
246  72  1.25E00  2.20E01 
249  72  6.65E00  2.85E01 
252  72  5.00E00  3.06E01 
255  72  1.25E00  5.00E00 
258  72  l.OOEOO  5.00E00 
261  72  1.25E00  5.00E00 
264  72  1.25 EOO  5.00E00 
267  72  6.61  E01  3.00E00 
270  72  1.66E01  5.00E00 
273  72  1.67E02  5.00E00 
276  72  l.OOEOO  5.00E00 
279  72  1.25E00  5.00E00 
282  72  1.25E00  5.00E00 
285  72  1.92E02  5.O0EOO 
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286  U  1.25E00  5.00E00 
289  U  7.64E00  5.00E00 
292  Tl 1.25E00  5.00E00 
295  Tl  6.60E00  5.14E00 
298  Tl 1.33E01  3.29E00 
301  Tl  1.25E00  5.00E00 
304  Tl 2.75E01  3.00E00 
307  Tl  3.24E01  3.00E00 
310  7e4.68E01  3.04E00 
313  71 8.25E00  4.00E00 
316  Tl  1.64E01  4.06E00 
319  n  1.25E00  5.00E00 
322  Tl  3.21  E01  6.30E00 
325  Tl  1.53E01  2.00E01 
328  Tl 4.05E02  4.91  EOO 
331  Tl  6.02E00  6.97E00 
334  Tl  1.52E06  3.97E00 
337  Tl  1.25E00  5.00E00 
340  Tl  3.00E01  5.00E00 
343  Tl  7.91  EOO  7.88E00 
346  Tl  8.85E00  6.64E00 
349  T  4.05E12  5.00E00 
352  T  7.41  El  7  1.72E01 
355  T  2.17E12  5.00E00 
358  T  2.05E12  5.00E0O 
361  T  1.05E14  8.52E00 
364  T  7.04E12  5.53E01 
367  Tl  1.25E00  5.00E00 
370  Tl 3.53E01  8.93E00 
373  Tl 1.50E01  3.00E00 
376  Tl 5.14E01  4.34E00 
379  Tl  1.36E01  5.00E00 
382  Tl  2.38E02  3.00E00 
385  Tl 2.99E01  5.00E00 


287  T  5.40E14  5.59E00 
290  Tl 1.25E00  2.64E01 
293  Tl  1.25E00  5.00E00 
296  Tl  1.25E00  5.00E00 
299  Tl  1.25E00  5.00E00 
302  Tl  1.55E01  5.00E00 
305  Tl  1.25 EOO  7.73 E01 
308  Tl  1.82E01  3.00E00 
311  Tl  1.20E01  5.00E00 
314  Tl 8.99E00  5.00E00 
317  Tl  1.68E01  5.00E00 
320  Tl  2.82E01  5.00E00 
323  ft7.80E01  3.16E00 
326  n  4.22E01  3.00E00 
329  Tl  1.25E00  1.25E02 
332  Tl  5.60E00  5.O0E00 
335  TZ  5.15E00  5.O0E00 
338  'll  1.25E00  5.00E00 
341  Tl  2.27E01  7.69E00 
344  n  9.99E00  4.90E00 
347  T  7.28E13  3.19E00 
350  T  1.30E10  5.00E00 
353  T  5.91  El 2  5.00E00 
356  T  1.75E12  5.00E00 
359  T  6.23E14  7.36E00 
362  T  5.60E14  6.23E00 
365  n  1.35E01  5.00E00 
368  U  1.25E00  5.00E00 
371  Tl  8.42E00  5.00E00 
374  U  1.02E01  5.00E00 
377  n  1.06E01  4.46E00 
380  n  1.25E00  5.00E00 
383  11 1.60E01  5.00E00 
386  U  1.25E00  5.00E00 


288  T  2.43E15  4.89E00 
291  Tl  4.24E02  4.30E00 
294  U  1.25E00  5.00E00 
297  n  1.20E01  1.21  E01 
300  U  1 .25 EOO  5.00E00 
303  11 1.25 EOO  6.91E01 
306  TZ  5.40E00  6.90E00 
309  U  5.84E00  1.33E01 
312  U  1.25E00  5.00E00 
315  7e3.76E01  4.10E00 
318  K  1.07E02  3.00E00 
321  U  2.37E01  3.50E00 
324  n  3.08E01  5.00E00 
327  'll  2.62E02  4.15E00 
330  Tl  1.25E00  5.00E00 
333  n  1.25 EOO  5.00E00 
336  n  7.18E00  5.00E00 
339  n  1.25 EOO  5.00E00 
342  Tl  7.61  EOO  5.00E00 
345  Tl  1.26E01  3.80E00 
348  T  8.63E11  5.00E00 
351  T4.99E11  5.00E00 
354  T4.83E12  5.00E00 
357  T  1.95E12  5.00E00 
360  T  7.09E14  4.56E00 
363  T  7.42E13  1.44E01 
366  7 1 1.25E00  5.00E00 
369  Tl  1.25E00  5.46E01 
372  Tl  2.06E01  3.00E00 
375  Tl  7.07E00  1.74E01 
378  Tl  1.25E00  5.00E00 
381  Tl  7.69E00  9.56E00 
384  U  1.19E01  5.00E00 
387  Tl  4.16E01  3.00E00 
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388  Tt  1.92E02  3.00E00 
391  Tt  5.30E00  5.00E00 
394  Tt  2.37E02  3.00E00 
397  Tt  1.25E00  5.00E00 
400  Tt  4.54E02  1.91E01 
403  Tt  1.25E00  5.00E00 
406  T  2.64E14  5.01  E00 
409  Tt  1.25E00  5.00E00 
412  Tt  1.25E00  5.00E00 
415  711.25E00  5.00E00 
418  Tt  l.OOEOO  5.00E00 
421  Tt  1.25E00  5.00E00 
424  Tt  l.OOEOO  5.00E00 
427  Tt  1.25E00  5.00E00 
430  Tt  l.OOEOO  5.00E00 
433  Tt  1.25EOO  5.00E00 
436  Tt  1.25E00  5.00E00 
439  Tt  1.25E00  5.00E00 
442  Tt  l.OOEOO  5.00E00 
445  Tt  1.80E03  3.00E00 
448  Tt  8.83E00  5.00EOO 
451  Tt  1.25E00  5.00E00 
454  Tt  1.25E00  5.00E00 
457  Tt  1.25E00  5.00E00 
460  Tt  1.25E00  5.00E00 
463  Tt  4.40E01  3.00E00 
466  Tt  l.OOEOO  5.00E00 
469  Tt  1.50E07  1.60E01 
472  Tt  l.OOEOO  5.00E00 
475  Tt  5.76E09  6.60E00 
478  Tt  l.OOEOO  5.00E00 
481  It  2.63E11  9.14E00 
484  Tt  l.OOEOO  5.00E00 
487  Tt  1.82E12  1.48E01 


389  Tt  1.25E00  5.00E00 
392  Tt  1.25E00  5.00E00 
395  Tt  1.25E00  6.69E01 
398  Tt  1.25E00  3.86E01 
401  Tt  1.25E00  5.00E00 
404  Tt  1.25E00  5.00E00 
407  T  1.75E14  4.84E00 
410  It  1.25E00  5.00E00 
413  Tt  1.25E00  5.00E00 
416  Tt  l.OOEOO  5.00E00 
419  Tt  1.25E00  5.00E00 
422  Tt  l.OOEOO  5.00E00 
425  Tt  1.25E00  5.00E00 
428  Tt  1.25E00  5.00E00 
431  Tt  1.25E00  5.00E00 
434  Tt  1.25E00  5.00E00 
437  Tt  1.25E00  5.00E00 
440  Tt  l.OOEOO  5.00E00 
443  Tt  1.80E03  3.00E00 
446  Tt  1.25E00  5.00E00 
449  Tt  1.25E00  5.00E00 
452  Tt  1.25E00  5.00E00 
455  Tt  1.25E00  5.00E00 
458  Tt  1.25E00  5.00E00 
461  Tt  1.25E00  5.00E00 
464  Tt  l.OOEOO  5.00E00 
467  Tt  7.63E05  2.35E01 
470  Tt  l.OOEOO  5.00E00 
473  Tt  1.13E09  7.64E00 
476  Tt  l.OOEOO  5.00E00 
479  Tt  8.58E10  7.81E00 
482  Tt  l.OOEOO  5.00E00 
485  Tt  1.25E12  1.26E01 
488  Tt  l.OOEOO  5.00E00 


390  Tt  6.97E01  3.00E00 
393  Tt  1.34E01  4.98E00 
396  Tt  1.25E00  5.00E00 
399  Tt  6.32E01  3.00E00 
402  Tt  1.25E00  5.00E00 
405  T  1.23E13  5.00E00 
408  T  7.86E13  4.66E00 
411  Tt  1.25E00  5.00E00 
414  Tt  l.OOEOO  5.00E00 
417  Tt  1.25E00  5.00E00 
420  Tt  l.OOEOO  5.00E00 
423  Tt  1.25E00  5.00E00 
426  Tt  1.25E00  5.00E00 
429  Tt  1.25E00  5.00E00 
432  Tt  l.OOEOO  5.00E00 
435  Tt  1.25E00  5.00E00 
438  Tt  l.OOEOO  5.00E00 
441  Tt  1.25E00  5.00E00 
444  Tt  1.25E00  5.00E00 
447  Tt  1.25E00  5.00E00 
450  Tt  8.83E00  5.00E00 
453  Tt  1.25E00  5.00E00 
456  Tt  1.25E00  5.00E00 
459  Tt  1.25E00  5.00E00 
462  Tt  1.25E00  5.00E00 
465  Tt  6.63E03  3.00E00 
468  Tt  l.OOEOO  5.00E00 
471  Tt  1.62E08  1.05E01 
474  Tt  l.OOEOO  5.00E00 
477  Tt  2.39E10  6.83E00 
480  Tt  l.OOEOO  5.00E00 
483  7e6.52Ell  1.06E01 
486  Tt  l.OOEOO  5.00E00 
489  Tt  1.96E12  1.85E01 
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490  TZ  1.00E00  5.00E00 
493  TZ  8.06 Ell  3.28E01 
496  TZ  l.OOEOO  5.00E00 
499  TZ  5.77E09  1.92E02 
502  TZ  l.OOEOO  5.00E00 
505  TZ  5.56E03  3.67E00 
508  TZ  2.31  E02  4.16E00 
511  TZ  4.20E01  4.25E02 
514  TZ  4.07E02  6.23E01 
517  TZ  1.25E00  5.00E00 
520  7 Z  l.OOEOO  5.00E00 
523  TZ  1.25E00  5.00E00 
526  TZ  l.OOEOO  5.00E00 
529  TZ  1.25E00  5.00E00 
532  TZ  l.OOEOO  5.00E00 
535  TZ  1.71E10  3.00E00 
538  7 Z  4.14E10  3.00E00 
541  7 Z  5.12E10  3.00E00 
544  TZ  5.62E10  3.00E00 
547  7 Z  3.93E10  3.00E00 
550  7 Z  4.64E10  3.00E00 
553  72.4.95E10  3.00E00 
556  7 3.96E10  3.00E00 
559  TZ  4.40E10  3.00E00 
562  TZ  4.54E10  3.00E00 
565  TZ  4.21E10  3.00E00 
568  72.4.34E10  3.00E00 
571  TZ  4.30E10  3.00E00 
574  TZ  4.25E10  3.00E00 
577  TZ  4.44E10  3.00E00 
580  TZ  4.20E10  3.00E00 
583  TZ  4.04E10  3.00E00 
586  TZ  4.73E10  3.00E00 
589  TZ  4.22E10  3.00E00 


491  TZ  1.52E12  2.37E01 
494  TZ  l.OOEOO  5.00E00 
497  TZ  3.39E10  1.49E02 
500  TZ  l.OOEOO  5.00E00 
503  TZ  1.15E04  4.38E00 
506  7e6.27E02  1.29E01 
509  TZ  8.20E01  6.90E01 
512  TZ  8.11E01  2.92E01 
515  TZ  1.25E00  5.00E00 
518  TZ  l.OOEOO  5.00E00 
521  TZ  1.25E00  5.00E00 
524  TZ  l.OOEOO  5.00E00 
527  TZ  1.25E00  5.00E00 
530  TZ  l.OOEOO  5.00E00 
533  TZ  1.25E00  5.00E00 
536  7£2.82E10  3.00E00 
539  TZ  4.55E10  3.00E00 
542  TZ  5.32E10  3.00E00 
545  TZ  2.84E10  3.00E00 
548  TZ  4.24E10  3.00E00 
551  'K4.77E10  3.00E00 
554  TZ  5.02E10  3.00E00 
557  7£4.17E10  3.00E00 
560  TZ  4.46E10  3.00E00 
563  TZ  4.56E10  3.00E00 
566  TZ  4.29E10  3.00E00 
569  72.4.33E10  3.00E00 
572  TZ  4.29E10  3.00E00 
575  TZ  4.65E10  3.00E00 
578  TZ  4.35E10  3.00E00 
581  TZ  4.14E10  3.00E00 
584  TZ  3.99E10  3.00E00 
587  72.4.52E10  3.00E00 
590  72.4.11E10  3.00E00 


492  TZ  l.OOEOO  5.O0E00 
495  TZ  2.61  Ell  5.45E01 
498  TZ  l.OOEOO  5.00E00 
501  TZ  1.92E09  5.72E01 
504  TZ  1.05E03  1.33E01 
507  TZ  3.56E02  1.25E01 
510  TZ  1.36E02  3.00E00 
513  TZ  1.75E04  3.00E00 
516  TZ  l.OOEOO  5.00E00 
519  TZ  1.25E00  5.00EOO 
522  TZ  l.OOEOO  5.00E00 
525  TZ  1.25E00  5.00E00 
528  TZ  l.OOEOO  5.00E00 
531  TZ  1.25E00  5.OOEO0 
534  TZ  l.OOEOO  5.00E00 
537  TZ  3.58E10  3.00E00 
540  TZ  4.87E10  3.OOE00 
543  TZ  5.48E10  3.00E00 
546  TZ  3.49E10  3.00E00 
549  TZ  4.47E10  3.00E00 
552  7£4.87E10  3.00E00 
555  TZ  3.63E10  3.00E00 
558  TZ  4.30E10  3.00E00 
561  TZ  4.51  E10  3.00E00 
564  TZ  4.58E10  3.00E00 
567  TZ  4.33E10  3.00E00 
570  TZ  4.32E10  3.00E00 
573  TZ  4.27E10  3.00E00 
576  TZ  4.54E10  3.00E00 
579  TZ  4.27E10  3.00E00 
582  TZ  4.09E10  3.00E00 
585  TZ  4.99E10  3.00E00 
588  TZ  4.36E10  3.00E00 
591  7e4.01E10  3.00EQ0 
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592  n  3.92E10  3.00EO0 
595  'll  5.27E10  3.00E00 
598  11 4.36E10  3.00E00 
601  1Z 3.89E10  3.00E00 
604  1Z  3.61  El 0  3.00E00 
607  1Z  4.64E10  3.00E00 
610  H  3.95E10  3.00E00 
613  ‘72.3.56E10  3.00E00 
616  ft5.10E10  3.00EOO 
619  7£4.09E10  3.00E00 
622  H  3.57E10  3.00E00 
625  11 5.83E10  3.00E00 
628  7£4.34E10  3.00E00 
631  11 3.64E10  3.00E00 
634  1Z  3.24E10  3.00E00 
637  11 1.73E11  9.74E00 
640  1Z  2.63E11  9.56E00 
643  1Z  3.19E11  9.43E00 
646  H  1.70E11  9.76E00 
649  1Z  2.51  Ell  9.53E00 
652  It  3.00E11  9.36E00 
655  11 1.72E11  9.82E00 
658  H  2.45E11  9.54E00 
661  1Z  2.87E11  9.34E00 
664  1Z  3.14E11  9.23E00 
667  7£2.45E11  9.55E00 
670  H  2.80E11  9.33E00 
673  1Z  3.02E11  9.19E00 
676  U  2.50E11  9.59E00 
679  1Z  2.77E1 1  9.35E00 
682  1Z  2.94E11  9.19E00 
685  1Z  2.60E11  9.71  E00 
688  U  2.79E11  9.37E00 
691  1Z  2.90E11  9.19E00 


593  H  3.85E10  3.00E00 
596  11 4.88E10  3.00E00 
599  72.4.17E10  3.00E00 
602  1Z  3.79E10  3.00E00 
605  7£5.49E10  3.00E00 
608  72.4.36E10  3.00E00 
611  1Z  3.80E10  3.00E00 
614  7£3.47E10  3.00E00 
617  714.67E10  3.00E00 
620  11 3.89E10  3.00E00 
623  11 3.45E10  3.00E00 
626  1Z  5.18E10  3.00E00 
629  U  4.06E10  3.00E00 
632  11 3.49E10  3.00E00 
635  11 7.31  El 0  9.97E00 
638  7£2.09E11  9.65E00 
641  11 2.85E11  9.49E00 
644  1Z  3.34E11  9.37E00 
647  1Z  2.03E11  9.65E00 
650  U  2.70E11  9.46E00 
653  1Z  3.12E11  9.33E00 
656TC2.02E11  9.71  EOO 
659  1Z  2.62E11  9.44E00 
662  1Z  2.97E11  9.31  EOO 
665  1Z  2.07E11  9.78E00 
668  1Z  2.59E11  9.45E00 
671  1Z  2.88E11  9.29E00 
674  1Z  3.07E11  9.18E00 
677  1Z  2.60E11  9.53E00 
680  1Z  2.83E11  9.31E00 
683  1Z  2.98E11  9.16E00 
686  1Z  2.68E11  9.56E00 
689  1Z  2.83E11  9.31  EOO 
692  11 2.93E11  9.14E00 


594  11 3.78E10  3.00E00 
597  11 4.59E10  3.00E00 
600  11 4.02E10  3.00E00 
603  n  3.69E10  3.00E00 
606  11 5.00E10  3.00E00 
609  11 4.13E10  3.00E00 
612  11 3.67E10  3.00E00 
615  11 5.68E10  3.00E00 
618  7£4.35E10  3.00E00 
621  11 3.72E10  3.00E00 
624  11 3.34E10  3.00E00 
627  H  4.70E10  3.00E00 
630  H  3.83E10  3.00E00 
633  11 3.35E10  3.00E00 
636  H  1.29E11  9.83E00 
639  11 2.38E11  9.62E00 
642  11 3.03E11  9.46E00 
645  H  1.28E11  9.92E00 
648  11 2.29E11  9.60E00 
651  1Z  2.86E11  9.41  EOO 
654  1Z  3.23E11  9.28E00 
657  U  2.26E11  9.60E00 
660  11 2.75E11  9.40E00 
663  1Z  3.06E11  9.26E00 
666  1Z  2.28E11  9.66E00 
669  H  2.70E11  9.40E00 
672  U  2.95E11  9.25E00 
675  U  2.36E11  9.73E00 
678  1Z  2.69E11  9.44E00 
681  U  2.89E11  9.24E00 
684  1Z  3.02E11  9. 11  EOO 
687  U  2.74E11  9.46E00 
690  1Z  2.87E11  9.24E00 
693  1Z  2.95E11  9.12E00 
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694  ft 2.97E11  9.09E00 
697  ft  2.85E11  9.43E00 
700  ft  2.89E11  9.23E00 
703  ft  2.93E11  9.07E00 
706  ft  2.96E11  9.52E00 
709  ft  2.92E11  9.26E00 
712  ft  2.91E11  9.08E00 
715  ft  3.15E11  9.60E00 
718  ft  2.99E11  9.31E00 
721  ft  2.92E11  9.10E00 
724  ft  2.87E11  9.00E00 
727  ft  3.10E11  9.38E00 
730  ft  2.96E11  9.13E00 
733  ft  2.87E11  9.00E00 
736  ft  3.23E11  3.00E00 
739  ft  5.82E11  3.00E00 
742  ft  7.29E11  3.00E00 
745  ft  3.23E11  3.00E00 
748  ft  5.58E11  3.00E00 
751  ft  6.84E11  3.00E00 
754  ft  7.62E11  3.00E00 
757  ft  5.51  Ell  3.00E00 
760  ft  6.56E11  3.00E00 
763  ft  7.20 Ell  3.00E00 
766  ft  5.59E11  3.00E00 
769  ft  6.44E11  3.00E00 
772  ft  6.92E11  3.00E00 
775  ft  5.83E11  3.00E00 
778  ft  6.44E11  3.00E00 
781  ft  6.77E11  3.00E00 
784  ft  6.98E11  3.00E00 
787  ft  6.57E11  3.00E00 
790  ft  6.72E11  3.O0EOO 
793  ft  6.82E11  3.00E00 


695  ft  2.81E11  9.67E00 
698  7£2.87E11  9.34E00 
701  ft  2.91E11  9.15E00 
704  ft  2.93E11  9.06E00 
707  ft  2.95E11  9.40E00 
710  ft  2.92E11  9.18E00 
713  ft  2.90E11  9.06E00 
716  ft  3.08E11  9.49E00 
719  72.2.96E11  9.23E00 
722  ft  2.90E11  9.07E00 
725  ft  3.28E11  9.60E00 
728  ft  3.05E11  9.26E00 
731  7£2.92E11  9.10E00 
734  ft  2.85E11  8.96E00 
737  7£4.29E11  3.00E00 
740  ft  6.39E11  3.00E00 
743  ft  7.65E11  3.00E00 
746  7£4.22E11  3.00E00 
749  ft  6.08E11  3.00E00 
752  ft  7.13E11  3.00E00 
755  7£4.29E11  3.00E00 
758  ft  5.93E11  3.00E00 
761  ft  6.81E11  3.00E00 
764  ft  7.35E11  3.00E00 
767  ft  5.93E11  3.00E00 
770  ft  6.63E11  3.00E00 
773  ft  7.04E11  3.00E00 
776  ft  6.08E11  3.00E00 
779  7e6.57Ell  3.00E00 
782  ft  6.85E11  3.00E00 
785  ft  6.40E11  3.00E00 
788  ft  6.63E11  3.00E00 
791  ft  6.76E11  3.00E00 
794  ft  6.84E11  3.00E00 


696  ft  2.83E11  9.54E00 
699  ft  2.88E11  9.28E00 
702  ft  2.92E11  9.11E00 
705  ft  2.99E11  9.64E00 
708  ft  2.93E11  9.33E00 
711  ft  2.91  Ell  9.14E00 
714  ft  2.90E11  9.02E00 
717  ft  3.03E11  9.39E00 
720  ft  2.94E11  9.16E00 
723  ft  2.89E11  9.01E00 
726  ft  3.18E11  9.47E00 
729  ft  3.00E11  9.19E00 
732  ft  2.89E11  9.06E00 
735  ft  1.86E11  3.00E00 
738  ft  5.13E11  3.00E00 
741  ft  6.88E11  3.00E00 
744  ft  7.96E11  3-OOEOO 
747  ft  4.98E11  3.00E00 
750  ft  6.49E11  3.00EOO 
753  ft  7.39E11  3.00E00 
756  ft  4.98E11  3.00E00 
759  ft  6.28E11  3.00E00 
762  ft  7.02E11  3.00E00 
765  ft  5.14E11  3.00E00 
768  ft  6.21E11  3.00E00 
771  ft  6.79E11  3.00E00 
774  ft  7.15E11  3.00E00 
777  ft  6.28E11  3.00E00 
780  ft  6.68E11  3.00E00 
783  ft  6.92E11  3.00E00 
786  ft  6.49E11  3.00E00 
789  ft  6.68E11  3.00E00 
792  ft  6.79E11  3.00E00 
795  ft  6.88E11  3.00E00 
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796  TZ  6.84E11  3.00E00 
799  TZ  6.77E11  3.00E00 
802  TZ  6.74E11  3.00E00 
805  72.7.30E11  3.00E00 
808  72.6.93E11  3.00E00 
811  TZ  6.74E11  3.00E00 
814  TZ 6.63E11  3.00E00 
817  TZ  7.20E11  3.00E00 
820  72.6.82E11  3.00E00 
823  TZ  6.60E11  3.00E00 
826  TZ  7.62E11  3.00EOO 
829  TZ 6.99E11  3.00E00 
832  7£6.63E11  3.00E00 
835  TZ  6.77E11  5.00E00 
838  TZ  1.31  Ell  1.38E01 
841  TZ  1.96E11  1.38E01 
844  TZ  2.35E11  1.38E01 
847  TZ  1.29E11  1.38E01 
850  TZ  1.87E11  1.38E01 
853  TZ  2.20E11  1.38E01 
856  TZ  1.31E11  1.38E01 
859  TZ  1.83E11  1.38E01 
862  TZ  2. 11  Ell  1.38E01 
865  TZ  2.29 Ell  1.37E01 
868  TZ  1.83E11  1.38E01 
871  ft2.06Ell  1.37E01 
874  TC2.20E11  1.37E01 
877  TZ  1.88E11  1.37E01 
880  TZ  2.04E11  1.37E01 
883  TZ  2.14E11  1.37E01 
886  7ei.98Ell  1.37E01 
889  TZ  2.07E11  1.36E01 
892  V,  2.12E11  1.36E01 
895  TZ  2.15E11  1.37E01 


797  TZ  6.81E11  3.00E00 
800  TZ  6.76E11  3.00E00 
803  TZ  6.74E11  3.00E00 
806  7£7.14E11  3.00E00 
809  TZ  6.85E11  3.00E00 
812  TZ  6.70E11  3.00E00 
815  TZ  7.66E11  3.00E00 
818  TZ  7.05E11  3.00E00 
821  TZ  6.74E11  3.00E00 
824  TZ  6.55E11  3.00E00 
827  TZ  7.36E11  3.00E00 
830  TZ  6.85E11  3.00E00 
833  7£6.55E11  3.00E00 
836  11 5.63E10  1.38E01 
839  ft  1.57E11  1.38E01 
842  TZ  2.11E11  1.38E01 
845  TZ  2.45E1 1  1.38E01 
848  R,  1.52E11  1.38E01 
851  ft2.00Ell  1.38E01 
854  TZ  2.28E11  1.38E01 
857  ft  1.53E11  1.37E01 
860  R.  1.94E11  1.38E01 
863  'll  2.17E11  1.38E01 
866  ft  1.58E11  1.37E01 
869  %  1.92E11  1.37E01 
872  TZ  2.11E11  1.37E01 
875  ?e2.23Ell  1.37E01 
878  TZ  1.95E11  1.37E01 
881  TZ  2.08E11  1.37E01 
884  TZ  2.17E11  1.37E01 
887  TZ  2.01E11  1.37E01 
890  ^2.09E11  1.36E01 
893  TZ  2.13E11  1.37E01 
896  TZ  2.14E11  1.36E01 


798  TZ  6.79E11  3.00E00 
801  TZ  6.75E11  3.00E00 
804  TZ  6.73E11  3.00EOO 
807  TZ  7.02E11  3.00E00 
810  TZ  6.79E11  3.00E00 
813  TZ  6.66E11  3.00E00 
816  TZ  7.40E11  3.00E00 
819  TZ  6.92E11  3.00E00 
822  TZ  6.66E11  3.00E00 
825  TZ  7.97E11  3.00E00 
828  TZ  7.15E11  3.O0E00 
831  TZ  6.73E11  3.00E00 
834  TZ  6.47E11  3.00E00 
837  TZ  9.81E10  1.38E01 
840  TZ  1.78E11  1.38E01 
843  TZ  2.24E11  1.38E01 
846  TZ  9.84E10  1.38E01 
849  TZ  1.71  Ell  1.38E01 
852  TZ  2.11  Ell  1.38E01 
855  TZ  2.35E11  1.38E01 
858  TZ  1.69E11  1.38E01 
861  TZ  2.03E11  1.38E01 
864  TZ  2.23E11  1.38E01 
867  TZ  1.72E11  1.37E01 
870  TZ  2.00E11  1.37E01 
873  TZ  2.16E11  1.37E01 
876  TZ  1.80E11  1.37E01 
879  TZ  2.00E11  1.37E01 
882  TZ  2.11E11  1.37E01 
885  TZ  2.19E11  1.37E01 
888  TZ  2.04E11  1.37E01 
891  TZ  2.10E11  1.37E01 
894  TZ  2.14E11  1.37E01 
897  TZ  2.13E11  1.36E01 
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898  ft  2.13E11  1.36E01 
901  ft  2.12E11  1.36E01 
904  ft 2.12E11  1.36E01 
907  ft  2.23E11  1.36E01 
910  7e2.15Ell  1.36E01 
913  ft  2.12E11  1.35E01 
916  ft  2.39E11  1.36E01 
919  2.22E11  1.36E01 

922  ft  2.13E11  1.36E01 
925  ft 2.08E11  1.35E01 
928  ft 2.31  Ell  1.36E01 
931  ft  2.17E11  1.35E01 
934  ft  2.08E11  1.35E01 
937  7 Z  1.25E00  5.00E00 
940  ft  1.25E00  5.O0E0O 
943  ft  1.25E00  5.00E00 
946  ft  1.25E00  5.00E00 
949  ft  1.25E00  5.00E00 
952  ft  2.32E06  3.00E00 
955  ft 1.25E00  5.00E00 
958  ft  5.08E00  5.00E00 
961  7ei.25E00  5.00E00 
964  ft  1.25E00  5.00E00 
967  ft  6.35E01  3.69E00 
970  ft  7.O0E00  3.10E01 
973  ft  2.10E02  3.37E00 
976  ft  1.21E03  3.00E00 
979  ft  1.01E03  3.00E00 
982  ft  1.43E02  3.00E00 
985  ft  1.72E01  3.55E00 
988  ft  2.19E08  2.17E02 
991  ft  1.25E00  5.00E00 
994  ft  1.25E00  5.00E00 
997  ft  6.16E00  5.00E00 


899  ft  2.12E11  1.37E01 
902  ft  2.12E11  1.36E01 
905  ft  2.12E11  1.36E01 
908  ft  2.20E11  1.36E01 
911  ft  2.14E11  1.36E01 
914  ft  2.11E11  1.36E01 
917  ft  2.31E11  1.36E01 
920  ft  2.18E11  1.36E01 
923  7£2.11E11  1.35E01 
926  ft  2.49E11  1.36E01 
929  ft  2.25E11  1.36E01 
932  72.2.13E11  1.35E01 
935  ft  2.06E11  1.35E01 
938  ft  1.25E00  5.00E00 
941  ft  5.36E00  4.97E00 
944  ft  1.25E00  5.00E00 
947  ft  2.21  E01  3.00E00 
950  ft  7.19E00  7.04E00 
953  ft  4.90E09  3.00E00 
956  ft  1.25E00  5.00E00 
959  ft  1.25E00  5.00E00 
962  ft  1.25E00  5.00E00 
965  ft  1.25E00  5.00E00 
968  ft  4.39E01  6.27E00 
971  ft  4.52E01  5.21  EOO 
974  ft  1.51E02  3.00E00 
977  ft  4.19E02  3.00E00 
980  ft  2.50E02  3.00E00 
983  ft  2.56E02  3.00E00 
986  ft  1.17E03  6.56E00 
989  ft  5.74E00  5.00E00 
992  ft  1 .25 EOO  5.00E00 
995  ft  1.25E00  5.00E00 
998  ft  1.25E00  5.00E00 


900  ft  2.12E11  1.36E01 
903  ft  2.12E11  1.36E01 
906  ft  2.27E11  1.36E01 
909  ft  2.17E11  1.36E01 
912  ft  2.13E11  1.36E01 
915  ft  2.10E11  1.36E01 
918  72.2.26E11  1.36E01 
921  72.2.15E11  1.36E01 
924  ft  2.09E11  1.36E01 
927  ft  2.39E11  1.36E01 
930  72.2.21E11  1.35E01 
933  72.2.10E11  1.36E01 
936  ft  2.12E11  1.36E01 
939  ft  1.25E00  5.00E00 
942  ft  1.25E00  5.00E00 
945  ft  1.25E00  5.00E00 
948  ft  1.25E00  5.00E00 
951  ft  2.19E03  3.58E00 
954  ft  8.30E00  5.00E00 
957  ft  1.25E00  5.00E00 
960  ft  1.25E00  5.00E00 
963  ft  1.25E00  5.00EO0 
966  ft  3.36E01  3.00E00 
969  ft  1.25E00  5.00E00 
972  ft  4.92E02  3.00E00 
975  ft  1.25E00  5.00E00 
978  ft  1.87E04  3.21E00 
981  ft  1.25E00  5.00E00 
984  ft  6.52E00  5.00E00 
987  ft  1.26E06  3.00E00 
990  ft  1.25E00  5.00E00 
993  ft  5.71  EOO  5.00E00 
996  ft  1.25E00  5.00E00 
999  ft  6.70E00  5.00E00 
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1000  ft  1.25  EOO  5.00E00 
1003  ft  8.05E01  3.00E00 
1006  ft 5.29E01  3.00E00 
1009  ft  3.62E01  1.62E01 
1012  ft  6.76E01  2.61E01 
1015  ft 5.56E02  3.00E00 
1018  ft  5.78E02  3.00E00 
1021  ft  1.25E00  5.00E00 
1024  ft  1.66E01  4.69E00 
1027  ft  1.25E00  5.00E00 
1030  ft 1.25E00  5.00E00 
1033  ft  3.75E01  5.00E00 
1036  ft  1.25E00  2.40E01 
1039  ft  1.25E00  1.54E02 
1042  ft  1.25E00  5.00E00 
1045  ft  1.25E00  5.00E00 
1048  ft  1.25E00  5.00E00 
1051  ft  1.00E00  5.00E00 
1054  ft  1.25E00  5.00E00 
1057  ft  1.25E00  5.00E00 
1060  ft  1.25E00  5.00E00 
1063  ft  1.25E00  5.00E00 
1066  ft  l.OOEOO  5.00E00 
1069  ft  1.25EOO  5.00E00 
1072  ft  1.25E00  5.00E00 
1075  ft  l.OOEOO  5.00E00 
1078  ft  1.25E00  5.00E00 
1081  ft  1.25E00  5.00E00 
1084  7 1  1.25E00  5.00E00 
1087  ft  1.25E00  5.00E00 
1090  ft  l.OOEOO  5.00E00 
1093  ft  l.OOEOO  5.00E00 
1096  T  8.54E18  4.59E00 
1099  T  2.15E20  5.21  EOO 


1001  ft  2.94E01  5.00E00 
1004  ft 1.25E00  5.00E00 
1007  ft  2.04E02  3.13E00 
1010  ft  1.25 EOO  5.00E00 
1013  ft  2.57E04  3.00E00 
1016  ft  1.25E00  5.O0EOO 
1019  ftl.25E00  5.00E00 
1022  ft  1.25 EOO  5.00E00 
1025  'll  3.66 E01  4.74E00 
1028  ft  1.25E00  5.00E00 
1031  n  1.25E00  5.00E00 
1034  n  2.09E02  3.00E00 
1037  n  6.72E02  3.00E00 
1040  T1.43E13  3.88E02 
1043  %  1.25E00  5.00E00 
1046  Tl  1.25E00  5.00E00 
1049  11 1.25E00  5.00EOO 
1052  %  1.25E00  5.00E00 
1055  n  l.OOEOO  5.00E00 
1058  U  l.OOEOO  5.00E00 
1061  n  1.25 EOO  5.00E00 
1064  11  1.25E00  5.00 EOO 
1067  n  1.25E00  5.00E00 
1070  U  l.OOEOO  5.00EO0 
1073  n  1.25E00  5.00 EOO 
1076  Tl  1.25 EOO  5.00 EOO 
1079  n  1.25 EOO  5.00 EOO 
1082  n  l.OOEOO  5.00 EOO 
1085  n  l.OOEOO  5.00 EOO 
1088  U  1 .25 EOO  5.00EO0 
1091  n  l.OOEOO  5.00E00 
1094  n  1.25E00  5.00E00 
1097  T4.42E19  5.05E00 
1100  n  1.25E00  5. 00 EOO 
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1002  n  6.86E01  3.00E00 
1005  n  7.78E01  3.00E00 
1008  Tl  7.43E01  9.76E00 
1011  H  2.06E02  2.41E01 
1014  Tl  1.61E03  3.00E00 
1017  TZ  3.45E02  3.00E00 
1020  K  1.25E00  5.00E00 
1023  U  1.25E00  5.00E00 
1026  H  2.03E02  4.57E00 
1029  U  1.25E00  5.00E00 
1032  n  1.76E01  5.00E00 
1035  n  1.05E03  5.00E00 
1038  n  9.21E02  5.00E00 
1041  n  1.25E00  5.00E00 
1044  n  1.25E00  5.00E00 
1047  Tl  l.OOEOO  5.00E00 
1050  ft  1.25E00  5.00E00 
1053  ft  1.25E00  5.00E00 
1056  ft  l.OOEOO  5.00E00 
1059  ft  1.25E00  5.00E00 
1062  ft  1.25E00  5.00E00 
1065  ft  1.25E00  5.00E00 
1068  ft  1.25E00  5.00E00 
1071  ft  1.25E00  5.00E00 
1074  ft  l.OOEOO  5.00E00 
1077  ft  l.OOEOO  5.00E00 
1080  ft  1.25E00  5.00E00 
1083  ft  l.OOEOO  5.00E00 
1086  ft  1.25E00  5.00E00 
1089  ft  1.25E00  5.00E00 
1092  ft  1.25E00  5.00E00 
1095  T  3.75E18  6.41  EOO 
1098  T  2.29E19  3.59E00 
1101  ft  1.25E00  5.00E00 
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1102  72  1.25E00  5.00E00 
1105  T1.35E13  5.00E00 
1108  T1.03E13  5.00E00 
1111  72  8.59E00  5.00E00 
1114  72 6.11E01  5.00E00 
1117  72 1.10E02  5.00E00 
1120  T8.26E11  5.00E00 
1123  T  1.85E11  5.00E00 
1126  T  8.92E10  4.99E00 
1129  72  1.25E00  3.43E01 
1132  724.15E01  3.88E00 
1135  72 7.57E01  3.00E00 
1138  T2.32E26  5.01  E00 
1141  T  5.20E25  5.00E00 
1144  T  2.51  E25  5.00E00 
1147  72  1.25E00  5.00E00 
1150  72 1.48E01  4.07E02 
1153  TZ  1.03E02  8.47E00 
1156  72  2.78E02  1.39E01 
1159  72  1.14E02  3.00E00 
1162  72  1.23E03  3.00E00 
1165  7?.  3.71  E03  3.00E00 
1168  TZ  1.07E03  3.00E00 
1171  TZ 4.59 E02  3.00E00 
1174  7 2  3.46E02  3.00E00 
1177  72  2.06E02  3.00E00 
1180  71  l.OOEOO  5.00E00 
1183  72 3.61  E02  3.00E00 
1186  71  1.90E04  3.00E00 
1189  7 2  1.72E02  3.00E00 
1192  7?.  2.83 E03  3.44E00 
1195  7?.  6.41  E02  6.96E00 
1198  TZ  1.48E03  4.23E00 
1201  71 7.25E02  3.00E00 


1103  TZ  1.25E00  5.00E00 
1106  T2.26E12  5.00E00 
1109  T  7.59E14  4.78E00 
1112  72. 2.06E01  5.00E00 
1115  72. 8.00E01  3.00E00 
1118  T  4.59E10  5.00E00 
1121  T  3.17E11  5.00E00 
1124  T  5.68E10  4.99E00 
1127  T  1.65E10  5.07E00 
1130  71 1.20E01  5.97E00 
1133  72. 5.54E01  3.40E00 
1136  T  1.29E25  5.01  EOO 
1139  T  8.92E25  5.O0E00 
1142  T  1.60E25  5.00E00 
1145  T4.64E24  5.00E00 
1148  TZ  2.00E00  5.00E00 
1151  TZ  1.25E00  5.00E00 
1154  TZ  1.14E03  5.34E00 
1157  TZ  7.19E00  1.03E02 
1160  TZ  1.23E03  3.00E00 
1163  TZ  1.97E02  3.00E00 
1166  TZ  1.65E03  3.00E00 
1169  72. 2.78E03  3.00E00 
1172  TZ  6.61E02  3.00E00 
1175  TZ  1.93E02  3.00E00 
1178  T  1.92E09  4.97E00 
1181  TZ  1.25E00  5.00E00 
1184  72.  2.00E04  3.00E00 
1187  TZ  1.18E03  3.00E00 
1190  72  7.19E03  3.00E00 
1193  TZ  1.43E02  4.15E00 
1196  72-  1.98E03  4.98E00 
1199  72-  6.31E03  4.03E00 
1202  72. 1.58E03  3.00E00 


1104  TZ  6.92E00  5.00E00 
1107  T2.06E12  5.00E00 
1110  TZ  2.32E01  3.41E01 
1113  TZ  4.12E01  5.00E00 
1116  72  9.65E01  3.00E00 
1119  T  1.71E11  5.01E00 
1122  T  1.41E11  4.99E00 
1125  T1.65E11  5.01E00 
1128  72. 1.68E02  3.00E00 
1131  72-  2.69E01  4.66E00 
1134  TZ  6.66E01  3.09E00 
1137  T4.82E25  5.00E00 
1140  T3.96E25  5.0OE00 
1143  T4.64E25  5.01EOO 
1146  TZ  l.OOEOO  5.00E00 
1149  TZ  2.27E02  3.00E00 
1152  TZ  1.82E03  7.09E00 
1155  TZ  1.42E01  6.65E00 
1158  72. 2.05E03  3.OOE0O 
1161  TZ  8.31E02  3.00E00 
1164  72  8.73E02  3.00E00 
1167  72  4.83E02  3.00E00 
1170  TZ  1.09E03  3.00E00 
1173  72  9.09E02  3.00E00 
1176  TZ  2.05E02  3.00E00 
1179  T  1.34E09  4.98E00 
1182  TZ  l.OOEOO  5.00E00 
1185  72  4.34E01  3.00E00 
1188  TZ  1.16E04  3.00E00 
1191  TZ  1.47E03  3.00E00 
1194  72  1.26E03  5.19E00 
1197  72  3.25E02  4.53E00 
1200  72  2.52E03  3.00E00 
1203  72  4.08E03  3.00E00 
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1204  n  1.30E03  3.72E00 
1207  U  1.03E03  4.02E00 
1210  U  1.41E02  6.85E00 
1213  T  3.76E23  5.01E00 
1216  T  3.32E12  5.00E00 
1219  T  4.09E09  4.99E00 
1222  T  9.35E26  5.00E00 
1225  T  1.15E24  5.01E00 
1228  Tl  1.25E00  5.00E00 
1231  %  5.63E00  5.00E00 
1234  U  1.25E00  5.00E00 
1237  T  1.76E11  5.00E00 
1240  n  1.25E00  5.00E00 
1243  n  7.50E00  5.00E00 
1246  T  1.19E08  5.00E00 
1249  T  1.52E11  5.00E00 
1252  1Z  6.65E04  3.00E00 
1255  1Z  1.36E01  1.38E01 
1258  n  1.25E00  2.45E01 
1261  U  8.61E00  1.27E01 
1264  n  8.97E02  3.00E00 
1267  1Z  4.11E03  3.28E00 
1270  'll  5.27E02  7.81E00 
1273  Tl  1.54E03  1.15E01 
1276  n  3.10E02  3.00E00 
1279  K  1.82E01  5.00E00 
1282  n  7.67E01  5.00E00 
1285  T  3.61  E14  5.85E00 
1288  T  1.33E13  2.66E01 
1291  T7.94E12  5.25E01 
1294  T  1.89E04  3.73E05 
1297  T  2.03 E04  1.10E05 
1300  T  7.19E02  5.00E00 
1303  T  5.98E09  1.59E06 


1205  n  5.41E02  3.81E00 
1208  ft  2.51  E02  6.54E00 
1211  1Z  1.41E02  6.86E00 
1214  T  1.97E09  5.01  E00 
1217  T  1.17E10  4.99E00 
1220  T4.15E09  4.99E00 
1223  T  3.29E24  5.00E00 
1226  T  1.17E24  4.99E00 
1229  U  1.37E01  5.00E00 
1232  n  2.15E01  5.00E00 
1235  T  3.07E13  5.00E00 
1238  T  8.85E12  5.00E00 
1241  U  1.79E01  5.00E00 
1244  Tl  1.29E01  5.00EOO 
1247  T  3.08E13  5.00E00 
1250  T  8.77E12  5.00E00 
1253  n  4.59E06  3.00E00 
1256  U  1.78E02  5.83E00 
1259  U  3.72E01  7.59E00 
1262  1Z  3.00E01  1.81E01 
1265  n  1.47E03  3.00E00 
1268  Tl  8.35 E01  5.35E00 
1271  11 1.10E02  1.71E01 
1274  11 1.36E02  5.00E00 
1277  U  3.53E02  3.00E00 
1280  11 2.99 E01  5.00E00 
1283  U  2.33E02  3.00E00 
1286  T1.04E15  3.60E00 
1289  T1.17E14  3.00E00 
1292  T  2.88E14  4.44E00 
1295  T  2.03E04  1.10E05 
1298  T  2.16E04  8.54E04 
1301  T  7.19E02  5.OOEO0 
1304  T  8.07E09  3.05E05 


1206  11 7.66E02  3.90E00 
1209  TZ  1.36E02  6.71E00 
1212  T  5.41E23  5.00E00 
1215  T  5.54E23  5.00E00 
1218  T9.03E09  5.01E00 
1221  T  1.70E06  5.00E00 
1224  T  2.54E24  5.00E00 
1227  T  4.79E20  5.00E00 
1230  11 5.31E00  5.00E00 
1233  11  1.20E01  5.00E00 
1236  T  8.35E09  5.00E00 
1239  T3.49E11  6.30E02 
1242  U  5.31E00  5.00E00 
1245  1Z  1.15E01  7.81E00 
1248  T  1.19E10  5.00E00 
1251  T5.59E11  5.00E00 
1254  11 1.05E01  1.05E01 
1257  11 2.13E02  6.45E00 
1260  n  2.66E01  9.02E00 
1263  TZ  6.66E01  1.54E01 
1266  1Z  3.11E03  3.17E00 
1269  11 8.92E02  6.02E00 
1272  1Z  8.95E02  1.10E01 
1275  U  1.94E02  5.00E00 
1278  U  4.54E01  5.00E00 
1281  11 8.35E01  5.00E00 
1284  T  2.41E14  5.97E00 
1287  T  1.27E15  3.13E00 
1290  T  2.41E14  4.03E00 
1293  T  8.37E14  4.50E00 
1296  T  2.1 6E04  8.54E04 
1299  T  1.15E03  5.00E00 
1302  T  2.40E10  6.96E05 
1305  T5.98E09  1.59E06 
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1306  T  8.07E09  3.05E05 
1309  T  1.43E17  4.99E00 
1312  T  2.16E04  7.53E04 
1315  T  2.15E02  5.O0EOO 
1318  T  1.92E04  8.70E11 
1321  T  2.04E04  4.65E11 
1324  T  7.49E09  5.00E00 


1307  T  2.64E17  4.99E00 
1310  T  1.89E04  4.01E05 
1313  T2.03E04  1.23E05 
1316  T2.10E02  5.00E00 
1319  T  2.04E04  4.65E11 
1322  T2.17E04  1.14E11 
1325  T  7.49E09  5.00E00 


1308  T  1.43E17  4.99E00 
1311  T2.03E04  1.23E05 
1314  T  2.16E04  7.53E04 
1317  T  2.10E02  5.00E00 
1320  T  2.17E04  1.14E11 
1323  T  5.32E09  5.00E00 
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